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SYMBOLS AND ABBREVIATIONS USED 


A 

a 

b 

bl 

c 

c* 

‘^eff 

F 

GEO 

So 

Hz 

h 

1 

I 

LEO 

M 

ID 

ID . 

,mol 

m 

^2 

Nu 

n 

ODE 

ODF 

ODK 

Pr 

P • 

R 

3^ 

Re 

St 


Surface area 
velocity of sound 
characteristic width 
boundary layer 
thrust coefficient 
velocity 

characteristic velocity 
effective discharge velocity 
thrust stress coefficient 
thrust; force 

geostationary orbit (2A hour orbit) 

acceleration due to gravity 

vibrations per second 

enthalpy 

inclination 

impulse 

low Earth orbit 
Mach number 

instantaneous mass of the gas storage container 
molecular weight 
Mass throughput 

molecular nitrogen (here purified) 

Nusselt number 

Direction of the normal (perpendicular) to a surface 
one-dimensional shifting equilibrium 
one-dimensional frozen equilibrium 
one-dimensional kinetic flow 
Prantl number 
pressure 

Special gas constant 
general gas constant 
Reynolds numerator 
Stanton number 
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s 

T 

TD 

TW 

t 

V 

X 

y 

z 


= coordinate tangential to the wall; entropy 
= temperature 
= two-dimensional 
= propulsive plant 
= time 
= velocity 

= nozzle axis coordinate 
= velocity distribution coordinate 
= altitude coordinate 


Subscripts 


acc 

ax 

c 

ch 

D 

1- d 

2- d 
del 
e 

F 

i 

inj 

lam 

n/v 

r 

sp 

t 

tan 

th 

turb 

ver 

vis 


= acceleration 
= axial 

= combustion chamber 

= constriction coefficient; throat throughput coefficient 
= characteristic 
= drag 

= one-dimensional 
= two-dimensional 
= delivered 

= end cross-section; end condition 
= thrust 

= initial condition 
= injector 
= laminar 
= inviscid 
= radial 
= specific 
= throat 
= tangential 
= theoretical 
= turbulent 
= vertical 
= viscous 
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W = balance 

® = environment condition (vacuum chamber or laboratory) 


Greek letters 


a 

a 

e 

“i 

V 
6 

6 * 

e 

n 

G 

\ 

X 

V 

V 
p 
T 

X 

vp 

ft 

0 ) 
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= heat transfer coefficient 
= exit cone angle (divergent portion) 

= exit cone angle (convergent portion) 
= isentropy exponent 
= boundary layer thickness 
= displacement density 
= energy density 
= surface ratio; eccentricity 
= viscosity 
r impulse loss density 
= thermal conductivity 
= geographical longitude 
= micro ( 1 0”^ ) 

= cinematic viscosity 
= density 
= thrust stress 
= compressibility 
= discharge coefficient 
= rising node 
= rotation rate 
= differential operator 
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1 INTRODUCTION 


/ 1 « 


At the present time, the areas of application of communications 
satellites include national and international telephone connec- 
tions, long-distance data transmission, radio and television 
signals, navigation for ships and aircraft, meteorological 
monitoring, satellite tracking and uses in development tasks. 

The optimum orbit is determined by the task posed. Synchronous 
satellites occupy an exceptional position, since they travel in 
their orbit at the same angular velocity as Earth (geostationary 
orbits). If such an orbit occupies the equatorial plane, then 
the subsatellite point remains always over the same location 
on the Earth’s surface: synchronous satellite. 

But conditions along a satellite orbit are not as ideal as one 
could hope. Reality is far removed from the assumption of a 
dynamic two-body system that can in addition be represented by 
point masses. Hence a satellite suffers a series of distortions 
such as those caused by 

- atmospheric resistance 

- anisotropies in Earth's gravitational potential 

- gravitational effects of other masses, especially the 
sun and the moon 

- radiation pressure due to direct solar irradiation. 

Earth radiation and albedo, as well as the satellite's 
own radiation 

- errors in the application of the control impulses 

- discharge of refuse, degassing, uncontrolled fuel re- 
leases 

- feedback coupling with attitude control systems 

- interactions with the environment (cosmic radiation, /2 
the Earth's magnetic field, etc.) 

- relativistic effects. 

* Numbers in the right margin indicate foreign pagination 
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Under the effect of these disturbances the satellite is not 
going to maintain its anticipated orbit. In the case of geo- 
stationary satellites this means that the subsatellite point 
no longer remains stable in one location, but wanders in an un- 
controlled manner. 

Among the disturbing effects listed above, three have proven to 
be predominant, for geostationary satellites: 

1) Anisotropy of the Earth’s potential 

2) Solar radiation pressure 

3) Gravitational effects of the sun and the moon 

Starting from the orbit distortions caused by these three in- 
fluences, the velocity adjustment can be calculated that is neces- 
sary to return to the initial position, once the satellite has 
achieved a known deviation from that position [ 13,3A,39,A2-AA ,66, 
79,81 ,82, 85-87,93,98, 110,129,131 ,152,1 6A, 166, 173-175) . 

The Tables below provide a rough overview on the order of mag- 
nitude of the gravitational distortion due to other bodies, on 
the geostationary orbit, and of that due to solar radiation 
pressure : 
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TABLE 1 


TABLE 2 


GROSS ACCELERATIONS DUE TO GRAVITATIONAL FORCES, 

OF BODIES IN THE PLANETARY SYSTEM, ON A SYNCHRONOUS 
SATELLITE 


Acceleration (m/s^) 


Sun 

5.931 

X 

10 

Mercury 

2.546 

X 

10 

Venus 

1.898 

X 

10 

Earth 

2.244 

X 

10 

Earth’s moon 

3.320 

X 

10 

Mars 

6.951 

X 

10 

Jupiter 

3.201 

X 

10 

Saturn 

2.324 

X 

10 

Uranus 

7.816 

X 

10 

Neptune 

3.708 

X 

10 

Pluto 

9.626 

X 

10 


-3 

-9 

-7 

-1 

-5 

-9 

-7 

-8 

-10 

-10 

-12 


COMPARISON OF THE PERTURBING ACCELERATIONS DUE 
TO SOLAR RADIATION PRESSURE AND DUE TO THE 
GRAVITATION OF THE SUN AND MOON 


Disturbance 

Solar gravitation 
Grav. due to moon 
Solar wind 
A/m = 2.05 m^/kg 
4.10 
6.40 
20.50 


Acceleration (m/s^) 


6.40 

X 

lO'^ 

2.44 

X 

10-^ 

6.25 

X 

10"^ 

1.25 

X 


1.95 

X 

10"^ 

6.25 

X 

10"^ 
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The distortions sDentioned are calculated with the aid of the /k 

variation in orbit elements, even though at point 1 the validity 
of the numerically determined results is limited, since the 
higher terms describing the Earth's potential field are affected 
by large uncertainties. The Earth's potential anisotropy and the 
sun's radiation pressure essentially cause distortions only in 
the orbital plane, while the influence on changes in inclination 
in relation to point 3, above, can be neglected. 

The consequence of equatorial ellipticity is a periodic motion 
around one of two stable points of the orbit, which are found on 
an extension of the equatorial ellipse's smaller semi-axis in 
both directions, at approximately 57° E and 123° W of geograph- 
ical longitude. The largest distortions occur at a satellite 
position A5° away from, such a stable point, attaining a radial 
deviation of ±Al.7km and an overall E-W deviation of 90°. The 
distortion's period is approximately 580 days. Due to the Earth's 
flattening, the radius of a synchronous orbit is increased 
marginally, by approximately 0.5 km. However, this is not relevant 
to orbit control. 

Orbit control becomes necessary when the satellite has reached 
a certain angular distance from its desired position. The return 
is accomplished by means of a classical Hohmann transition. The 
necessary velocity demand is independent of the number of control 
maneuvers but depends on the satellite's initial position with 
respect to the closest stable point. If this distance is A5°, 
then the velocity demand attains a maximum of 

AVpo^ = 5.17 m/sec/yr. 

How large the distortion due to solar wind results for the 
satellite depends on the ratio of its cross-sectional area A 
and its mass m. 

For A/m = 10"3 m‘^/kg, the maximum angular deviation from an ar- /5 
bitrary initial point is AA = 0.25° and a maximum radial deviation 
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is r = 17.7 km. If an angular deviation of 0.1® is allowed, as 
is customary with current communications satellites, then a 
yearly velocity demand of 


= 3.75 m/sec/yr 


must be satisfied. 

The gravitational forces of the sun and the moon (point 3) 
cause a distortion normal to the orbit's direction that becomes 
manifest as a change in inclination. This distortion has a period 
of p = 53.6 yr and a maximum deviation of degrees. The 

initial change of inclination occurs at Ai = 0.865®/yr. Applying 
a simple impulse, after a certain time and angular change the 
satellite is returned to its old position. During a five-year 
mission, with a permissible inclination change of 0.25®, the 
velocity demand is 


Av = 48.55 m/sec/yr. 
grav 

By applying the "phase angle technique" this value can be reduced 
to 


Av = 45.85 m/sec/yr. 
grav 

Normally orbital control is accomplished by means of chemical 
propulsion systems. This does not apply to the attitude control 
of geostationary satellites, which generally are equipped with 
a combined system of mechanical devices (gyroscopes) and reaction 
equipment. For this reason it is not possible to provide a 
precise value for the velocity demand for attitude control. At /6 
a permissible deviation of 0.1® a realistic value of 

Av_. . = 1 m/sec/yr 

O C Wr 

may be assumed. 
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Thus, if a satellite is on a weakly elliptical orbit of 
inclination i = 0, on a synchronous Earth orbit, with an 
allowed longitudinal deviation of 0.1® and a 0.25® allowed 
lateral deflection, then a yearly velocity demand of at least 

AVtot = 55.77 m/sec/yr 

must be provided for, assuming the orbital distortions are 
additive, in regard to the Av due the sun, moon and Earth. 

The requirement for high precision in maintaining a satellite's 
orbital position leads to the installation of orbit and attitude 
control systems. 

Orbit control compensates for center-of-mass distortions, while 
the function of attitude control is to insure that the satellite 
retains its alignment with a certain target point. A system that 
satisfies both functions is highly desirable, but this goal has 
not been achieved, as yet. 

In recent years a number of chemical and electrical [21,92,151] 
propulsion systems have been developed (as well as gyroscope 
systems for attitude control) and successfully installed in 
operational satellites. Many others are still in the planning 
stage. Consideration was also given to using the solar wind or 
the light pressure for orbit and attitude control. The extensive 
activities precisely in the filed of control are not surprising, /7 
since due to the high cost of space travels any possibility for 
savings must be taken advantage of. 

One look at the basic equation of rocketry, Av = c ^-*ln (m /m ) 

et I o e 

makes it very clear that for a given velocity demand and mini- 
mization of the fuel mass (m^ - the highest possible exhaust 

velocity must be adopted. However, the consequent demand for the 
use of stable two-component systems of intermediate energy meets 
a series of technological restrictions: the required control 
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precision of the satellite with regard to orbit and attitude 

leads to the application of a large number of correction 

maneuvers with individual burning durations (typically 10,000 

per year). Simultaneously, in order not to reduce the operating 

time per propulsion impulse to unrealistically brief periods 

of time, the propulsion unit’s thrust must be correspondingly 

low. While for N/S corrections two-component systems can still 

be used (for instance, MBB propulsion unit at 10 Newton; = 

290 sec), there are limits on miniaturization in the range of 

0.5 to 2 N. Today, ^ 2 ^^ single-component fuel systems are used 

here (without heating, I = 230 sec; with electric heating: 

sp 

I = 295 sec). Even smaller thrusts (50 mN ) , as required by 
sp 

precision attitude control, are accomplished by means of cold 
gas propulsion units (for instance, He of GN^ [gaseous N^l), 
since even the catalytic hydrazine devices can not be arbitrarily 
miniaturized. The low specific impulse of the engines (for nitro- 
gen, typically 70 sec) forces an optimization of them with respect 
to combustion chamber pressure, expansion ratio, number of oper- 
ating periods and overall consumption. 

It is the purpose of this investigation to provide a contribu- 
tion in this regard. 

During the execution of the PAEHT study [56-63], the lack of 
computer programs adequate to precise performance prediction of 
chemical rocket propulsion units became very clear, in Europe 
[2,6,10,18,32,36,37,^8,76,77,83,85,86,88,109,116,120,121,130, 

132, 137 , 138, 148, 149, 156, 162] . The programs developed in the US /8 
under the superordinate concept JANNAF Procedures (ODE, ODK, 

TDK, BLIMP) could be obtained but were only partially usable (not 
suf f iciently ) in Europe [16,22]. 

For these reasons the current study was initiated at the con- 
clusion of the studies performed for the PAEHT Project by the 
Chair for Space Technology of the Technical University Munich. 
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The first main portion of the task consisted in the development 
of a simplified, explicit procedure for performance prediction, 
in the foreground, and the comparison with available theoretical 
results from the procedures mentioned above. On the basis of 
these results, the final objective was a collection of easily 
manipulated estimation formulas, such as those required in the 
design of rocket propulsion plants. Naturally, this also requires 
the corresponding comparisons with measured values. However, 
during the course of the studies already mentioned, considerable 
gaps were found also in the area of reliable, available perfor- 
mance measurement values for the case of small propulsion units 
as they are particularly in the foreground, here. 

This gap is partially closed by means of the test program related 
to cold gas propulsion units, presented in Section 2. 
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2 THEORETICAL INVESTIGATIONS ON PERFORMANCE PREDICTION 


/9 


2.1 Prediction by extrapolation of data from implemented pro- 
pulsion units 

There are two methods, in principle, for the prediction of per- 
formance parameters (not only of propulsion units) of future 
systems: extrapolation, starting from data of actually implemented 
equipment, or by determining the physical laws and interactions 
that govern the individual components. 

The first of these requires extensive data collections for 
implemented systems, including parallel developments of different 
devices. Only thus can hasty judgements be prevented. The ad- 
vantage of this way is that it is not necessary to delve very 
deeply into the individual systems and the physical laws that 
govern them, provided it is assured that the groups or companies 
participating have or had available the necessary experience. As 
applied to the chemical propulsion units (more especially, of 
low thrust) currently at the focus of our interest, however, the 
following serious problems occur, which we shall mention here 
only in terms of keywords: 

- large numbers of different concepts 

- few parallel developments 

- To date, no large areas have been covered; rather, almost 
always special applications are placed in the foreground 
of a given development 

- large differences in the technical implementation 

- performance differences among participating groups 

- insufficient documentation of the results. 

For this reason we have preferred following the second way. 
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2.2 Performance prediction by determination of individual /TO 

losses 


The starting point for the determination of individual losses 
in a rocket propulsion unit is the zero-loss, one-£imensional 
•'shifting equilibrium" flow [ODE], whose performance values 
represent the upper theoretical limit. It is assumed that the 
gas composition corresponds in each case to the pressure and 
temperature conditions during expansion in the nozzle, and that 
the necessary reactions occur without delay time. The energy 
equation provides the isehtropic exhaust velocity: 


‘sp.ODE.vac So • V'2(h^ - h^) with 'c ' ‘p.e 


(I) 


However, in the real flow a series of performance-reducing 
factors occur, as a consequence of departure from the ideal 
condicions, I^p (see also Figures 1, below and 2.1) 


Injection zone 
(Mixing effects 
predominate) 


Displacement thickness 6* 
Impulse loss thickness 


Boundary layer 
thickness 



Fuel 


sp.del 


Injector 


'Mixing Evaporation effects range 
& evap. predominate; begin 
effects of boundary layer 
formation 


j expansion; , |Maptacion' 

I 2-dimension.' to environ- 
reaction ki- mental con- 
netics&bdry ditions 
layer effects 
predominate 


Figure 1 Processes in the combustion chamber and nozzle 
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^sp.del “ ^sp,ODE,v«c ’ ?er * \in * *^part * ^'Z-d.O-d) 


( 2 ) 


'bl 


^adapt 


'»h 


2.2.1 Energy conversion in the combustion chamber /II 

The first efficiency term on the right side is called energy 
conversion efficiency (energy release) of the injector-combus- 
tion chamber system. In the general case of the chemical multi- 
component rocket propulsion unit it is not possible to obtain 
the total, theoretically possible energy yield, starting from 
the chemicals injected into the combustion chamber. Because of 
the finite dwell time in the combustion chamber (v^ / 0), in- 
complete mixing and fuel preparation, as well as due to technical 
implementation restrictions, can not become equal to 1.0. 

However, implemented systems (SSME) today reach 0. 997 . Since here 
GN 2 occurs only as a single-component system and no reaction 
takes place, we can here set = 1.0, however. 

2.2.2 Kinetic efficiency 

During the gas expansion process within the nozzle flow, the 
thermodynamically limiting parameters - such as T, p, p, X, n - 
undergo significant changes and hence the combustion chamber 
conditions of the chemical condition ODE are no longer appli- 
cable. For instance, at high dissociate molecules tend to 
recombine, as T decreases. Similarly, charged atoms and molecules 
(ions) are converted into neutral particles, by charge exchange. 

In addition, the state of excitation in the combustion chamber 
can initiate vibration processes in molecules, given a sufficently 
high Ah. All of these phenomena have in common that if the pres- 
sure or the temperature of the working fluid decrease, given the 
availability of a sufficiently long period of time, recombination 
or respectively, a relaxation of the excited degrees of freedom 
would occur. However, if this time - which depends very critically 
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on the gas component’s characteristics, in each case - is not 
made possible by assuming an infinitely large expansion 
velocity, then we refer to a one-dimensional frozen equilibrium 
[ODF]. It represents a lower theoretical chem-thermodynamic /^2 
performance value, because the energy contained in the gas' 
dissociation, ionization and vibration is not available - in 
contrast to the ODE flow - for expansion work and hence, the 
required gas acceleration; rather, it is still contained in the 
gas stream beyond the nozzle, after leaving it, in the form of 
internal energy, without contributing positively to the thrust. 

In the "technical nozzle flow" the so-called ODK flow [one- 
dimensional kinetic equilibrium] always lies between ODE and 
ODF (see Figure 2.2). 


Hence in general the dwell time is sufficient to include a part 
(but of course not all) of the difference h(ODE-ODF) in the 
expansion work. The lower kinetic performance value can be 
expressed as 


I 


sp, ODF, vac 


^ODF PDF. va c 

9o 


(3) 


Usually, to facilitate individual calculations, the total 
kinetic efficiency can be broken down into the product 

^kin * \,dis ’ \,ion ’ ^k,vib 


At this point it must be recalled, in general, that in the area 
of chemical kinetics considerable gaps still exist in the sense 
of an applicable engineering approximation. For some very special 
cases (molecules or gases, respectively, at certain specific 
pressure and temperature conditions) there are indications, in 
the literature, on reaction and recombination constants, but the 
quantitative theoretical statements contradict each other [71] 
as well as the few available measurement results. As an example 
we shall consider only the simplest case of the biatomic gas N 2 
in the excited state. 
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According to Landau and Teller [53] the energy exchange between 
the vibrational and translational degrees of freedom can be 
described by means of the following linear differential equa- 
tion {based on the model of the harmonic oscillator) 

° _ ®v,Equi “ ®v 

ITT. T; (5) 

where e^ is the vibrational energy per unit mass of gas; 
e pnni value of e when the molecule's vibration and 

translation are at equilibrium; and is the relaxation time 
for the exchange of vibrational and translational energy. 

Assuming a Boltzmann distribution for the vibration energy steps, 
vibration temperatures can be calculated for them and then e^ 
will be determined by 


,, Rh 
* k - I) 


( 6 ) 


Hence can be derived using the local pressure p as parameter 
and compared to measured experimental data, for instance with 
the aid of Raman spectroscopy. The results of this comparison 
(the experiments were performed several times, by different 
authors) are a surprise until this day, and are still awaiting 
clarification. According to Landau-Teller the time would be 
higher by a factor of 70 than the measured relaxation time. 


A 


Sebacher hence foregoes a theoretically derived function and 
determined a compensating function adapted purely to empirical 
values : 

\ P - 3 • 10*’2 ^ (bar sec) - (7) 

which here is typically 10 sec (T is the absolute gas tempera- 
ture, here, NOT to' be confused with T^ above). 
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Figure 2.1 Qualitative course of the impulse, taking the various 
losses into consideration 

Comparison with the residence time of the gas in the propulsion 
unit 


,c Fx0,5 


typically io"^ sec 


( 8 ) 


shows that in practice the flow becomes "frozen” shortly behind 
the throat area, such that in the recombination equation 


I 


sp,0DF,vac 


- k ♦ (1 - k ) 

“ *sp.0DE,vdC 


( 9 ) 


can be set equal to 0, and we thus finally arrive at = 1.0 /15 

for the case here under consideration. 


This allows us to complete all remaining comparison calculations 
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E — 

Figure 2.2 Qualitative courses of ODE, ODK and ODF impulses, as veil as 
kinetic efficiency, as a function of the surface ratio [47] 

with ODF parameters, which is very advantageous because then 
subsequent performance calculations can be separated from the 
kinetics problem. It is certain that also in the future chemical 
kinetics will considerably complicate total performance calcula- 
tions for more complex gas combinations [33,91,112,125,138,157- 
160] 


2.2.3 Particle flow 


/16 


In the case of the monopropulsion unit, gas particle flow as it 
occurs in modern solid fuel propulsion units (high A1 proportion) 
, as well as condensed droplets, such as they occur in liquid 
fuel engines, do not take place. Nor was it possible to observe 
any wear in the galvanized nozzle wall (not even in the espe- 


cially endangered throat region). Hence, 


r)part“^*° shall be used. 
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2.2. h Multidimensional flow losses 


2. 2.^.1 Divergence losses 



Figure 3 Geometry model to determine divergence losses 


Due to radial divergence of the flowlines in the nozzle cross- 
section as well as the distribution of the velocity and pressure 
across the end cross-section surface area, there is an impulse 
loss (in relation to l-d flow). (3-d effects, such as vortex 
formations due to asymmetry of the propulsion unit geometry, in 
general can be disregarded as very small). In the case of the 
bell nozzle, these multidimensional flow processes are usually 
numerically simulated by means of characteristics procedures 
[1A4]. In order to estimate ID losses in conical nozzles, how- 
ever, the really rather voluminous computer programs are not 
necessary. Starting from a source flow (such as, for instance, 
the potential theory, [161]), the thrust fraction of the annular 
spherical surface element dA^ in Figure 3, above , [ 1 1 9 , 1 79 ] is 

dF « dm cos 0 ♦ (p^ - p^) dA^ cos G /17 

With dm » snd “ 2 uR* sin G d 0 and integrating once, 

Fjo • 2 TiR*r ♦ (p^ - pj 

•'o 


] 


cos G sin G d G 


( 11 ) 
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Assuming that p , v and p are constant over the spherical 

s s s 

surface A^, and with 


A ■ - 2 tiR* cos of* ■ 2 nR* (1 - cos a ) 
* 0 * 

(Ap ■ tiR* sin* a^) 


(12) 


equation (11) is rewritten as 



1 ♦ cos 


(- 


T 


) ^ ^'s''s ^ 


(13) 


With in “ Pj Aj the expression between brackets becomes 

Fj « m ♦ (p^ - p^) (14) 

F is the 1-dimensional thrust derived from the source flow 
s 

model. The divergence factor corrects the 1-dimensional 
thrust taking into account the divergence losses and depends /18 
only on the opening angle of the divergent portion of the 
nozzle : 

njjj • 0,5 (1 ♦ cos Og) 


and hence 


^TD " *^s 

With the aid of a Fortran program [1^^] available at the Chair 
for Space Technology it could be shown that the simplifying 
assumptions that led to the determined above for conical - 
nozzles, are justified. The differences between the nyo factor 
determined explicitly by means of the characteristics procedure 
and the simplified theory is less than 0.3% for a 15° conical 
nozzle [ 1A,51 ,5^,72,73,89,95,112,122-124, 135, 138, 14A,U7, 154,157) . 
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Conical nozzle 


Figure 4 Two-dimensional flow losses for bell and conical 
nozzles as a function of normalized nozzle length 


Figure 4, above, shows typical comparative results of the 
calculation for conical and bell nozzles. 

2. 2. A. 2 Mass throughput reduction /19 

Because of constriction and boundary layer effects the actual 
mass throughput is reduced, with respect to 

In the LRT characteristics program already mentioned, a TD throat 
isobar is determined, starting from a Sauer-Nicholson procedure, 
which actually serves as a starting point for the entire pro- 
cedure. It is not only possible to optimize the TD flow, thereby 
(including wall pressure distribution), but a TD-m-eff iciency 
can also be calculated, which usually lies between 0.993 and 
0.996 (it does not include friction or boundary layer influences, 
yet). Based on the work of Back et al . [3-5], for instance, in 

addition the dependence of n* on the throat’s radius of curvature 

m 

can be determined. The result is shown in Figure 5 (page 27), 
normalized with respect to r^. It can be seen that the lowest m 
loss occurs when the throat's transition radius r^ is equal to 
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Figure 5 Effect of the throat transition radius 
rj. on the TD losses 127,69] 


r^. The most substantial decrease in the ODE mass flow can be 
attributed to viscous distortion effects, which "clog" the 
nozzle throat area. 


Hall et al [53] proposed the following equations for discussion 


m. 


'O " 




'ODE 


r. 2 
(t ♦ 1)0 
r 



8 Y ♦ 21 r. 


754 Y* ♦ 1971 Y ♦ 2007 
55295D 



] 


(17) 


Sauer, Oswatitsch and RothsteinL 11 A , 133 , 1 34 ] had already tried 
a similar equation, which does not require the second and third 
terms in the square bracket. Kuluva and Hosack [89] propose a 
somewhat different representation, 


r + 0.05 r. 0,019 r ♦ 0.1 r. 0,21 , 0,5 t 

C. 7 S r^ > • <Ri> -WJ. (16) 


which in a second range, of 30 < Re < 10^ agrees well with the 
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. experimental data. Figure 6, below, provides an overview; it 



Throat Reynolds number 

Figure 6 Mass throughput coefficient c^ as a function of Reynolds number 
in the throat. Experimental data: r^/rt = 4.42; y = 1.4. Solid 
curve using equation (18) 

can be seen that even at Re = 2,300, = 0.90, which differs 

from other results reported. All of the statements found in the 
literature have in common that they were obtained with substan- 
tially larger nozzles (throat diameters) and a stronger influence 
of local heat flows. In the case of the small diameter here 
investigated, with a throat of 0.2 mm, and adiabatic conditions 
in the flow range along the walls, the characteristic throughput 
values depart from the equation above (see section 4.2). As was 
the case in the correction procedure mentioned in section 
2. 2. 5. 3 for the determination of u^, in the determination of Cj^ 
we can resort directly to the available boundary layer quantity 
6* and 0. Since here - in contrast to nearly all procedures 
reported in the literature - in the area of the throat 6* ^ 0 

and e 0 (see also section 2.2,8), the contraction (or throat / 2 ^ 

A* 

throughput) coefficient ^ - t«v can be calculated with 

° *t,n/v.0DF 

sufficient accuracy by numerical methods. In anticipation of 
sections 4.2 and 5, this way proves to be fully justified. 
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2.2.5 Boundary layer loss calculation 


Because of gas viscosity and surface roughness on the propulsion 
unit’s inner walls, boundary layers are formed, with the effect 
of reducing performance. (In large propulsion units these losses 
are merely a few (1-2) %, while they are substantially higher 
for smaller engines.) The first decision in the calculation of 
boundary layer problems requires an answer to the question of 
whether the flow is laminar or turbulent (see also [11,2^,26,29, 
30,38,A1 ,65, 9^,95,99,100, 103, 104, 107, 108, 112-1 U, 11 7, 118, 122-124, 
126,135,157,170,172,177,179)). 

2.2.5. 1 Flow condition 

Normally, laminar flow can be assumed for nozzles of small 
dimensions. Because of the high flow velocities, the large 
Reynolds numbers would permit turbulent flow and it hence must 
be shown whether the large acceleration causes a laminarization 
of the boundary layer, within the nozzle. According to Boldman 
et al. [15], the following criterion can be used; 


K 


acc 


ar * ^ 23^) > 2.88 • 10"® 


(19) 


It can be readily shown that in the case of small nozzles the 
acceleration parameter K _ is several-fold above the critical 

aG C 

value. Hence, all boundary layer calculations presented will be 
based on a laminar velocity profile. 


2. 2. 5-2 Impulse equation 

The starting point for boundary layer calculations is the impulse 
equation for axisymmetrical nozzle flow [1,7-9,96,97,136,146,167] 


do 


"f 

T 



"2 - M' ♦ 6*/G 

dM 1 dr*] 

|l . (jj) J -0 

.(1 K*)M 

’ 3z *?3lJ 


( 20 ) 


/22 


where 
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is the displacement thickness of the boundary layer (Figure 7): 

/23 



Kozzle axis 

Figure 7 Designation of main flow lines 



is the impulse loss thickness, and 





( 22 ) 


( 23 ) 


the dimensionless friction coefficient with respect to wall 
shearing stress, M the Mach number of the core flow, r the local 
nozzle radius and d/dz the differential with respect to the 
nozzle axis coordinate. 


In order to calculate the Mach number, we shall resort to an 
equation that can be solved only by iteration, for a given 
surface ratio e: 
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2. 2. 5. 2.1 Geometry determination 

For the boundary layer calculation to be sufficiently precise, 
it is necessary to simulate the exact nozzle contour, without 
edges or obstructions. In the EGVP program here developed, the 
task (R = f(z)) is solved by a subroutine in six sections, wit 
the aid of the main geometry points shown in Figure 8, below. 


Conical nozzle throat 


Figure 8 Nozzle geonetry 


1st section: ZO to ZC 
2nd section: ZC to B1 
3rd section: B1 to B2 


4th section: B2 to H8 
5th section; H8 to B3 
6th section: B3 to ZE 

In contrast to most authors, who start the boundary layer 
calculations at the nozzle's narrowest cross-section - which 
seems justified in the case of large propulsion units, because 
the severe contraction of the flow at the throat reduces the 
boundary layer thickness to practically 0, at the narrowest 
cross-section - here we begin at the combustion chamber head. 

It is shown that even for small throat diameters (< 1 mm) the 
boundary layer thickness appears to be negligible in the throat 
region (see Figure 9, below). Nevertheless, including the subsonic 



m 

3 

(0 

u 

0) 


M 

N 

o 

c 


cd 

V 

o 

*-3 


vail 


range in the boundary layer calculation of nozzles represents a 
physical expansion and numerical completion - 
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2. 2. 5. 2. 2 Shearing stress coefficient Cj. /26 

In general, for locally constant viscosity the shearing stress 
at the wall is 

\ ( 26 ) 

where is the gas viscosity at the wall, which is a function 
of the temperature: 

'’w • ^ (^w^ - f (M. Y, e) (27) 


The shearing stress coefficient is defined as 


‘f 



^ ^ 5y 


( 28 ) 


Figures 10 and 11, below, show the nomenclature, in further 



Figure 10 Linear velocity 
profile 



profile 


considerations, for a general and a linear velocity profile 
within the boundary layer. 

The gas density at the wall is derived from 


/27 
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1 ) 0 ) 


'c 


Y-1 


1 ♦ — g — 


Mi 


and 


(31) 


u = u = 0 we can determine p , the density at the wall, 
w w 


“ Pc 


y=T 


(32) 


where the subindices refer to: 
c = chamber conditions, 

g = edge of the boundary layer (to the core flow), 
w = wall . 


We must now take into account the temperature dependence of the 
viscosity, r\ = ri(T). Anticipating experimental results, we 
shall assume that no heat is removed through the wall. With 


T 


■ 1 


_ 


«9 


. - 


and 


(33) 


u ■ ■ 0; \ = 0 we obtain 


■ Tg and hence also ^ \ 


(34) 


With equations (29) and (32), the shearing stress coefficient Cj. 
becomes : 


3A 
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( 35 ) 





The temperature dependence of the GN^ viscosity is approximated 
by means of 


where 


n 



(36) 


Tg • 300 K 

^ kg/msec 

X • 0,71976 

A variety of expressions is used in the literature: 


"f 


0,5533 


. 0.571 


(Ziebland (177)) 

(37) 

(Bartz (9)) 

(38) 

(general; see Figure 

12;[61]) 
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Here Reg represents the Reynolds number related to the impulse 
loss thickness. 


c 




(39) 
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(SchUchtlng (136)} 


'f 





c. . irMi - M12§§ (Shapiro (150)) 

^ >/Ri; “ 



Figure 12 Effect of the shearing stress coefficient K] 
on the boundary layer losses 


2. 2. 5. 2. 3 Derivation of the boundary layer shape factor, 6*/0 

Using the usual integral equations of the boundary layer theory 
for the displacement and impulse loss thicknesses, 



where p and u are the density within the boundary layer, p 

O 

and u the corresponding values at its edge, and ^ is the 


/30 


36 



coordinate perpendicular to the wall, the ratio 6*/6 is 
sought, which is necessary for the integration of the impulse 
equation. 


By introducing the ideal gas equation and using equation (38), 
it is possible to rewrite equation (29) thus 


p 





1 


T ^ 

c '9 'c . /U / 

eg g 


(44) 


Using the linear velocity distribution on which equation (36) 
is based (a comparison with the exact distribution is shown in 
Figure 13, below), the integral equation can be rewritten: 



Boundary layer edge 


Figure 13 Comparison between the simplified linear and actual velocity 
distribution with respect to the course of density and 
velocity within the boundary layer 
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ly/g) 




and 


(45) 


37 


, where (^6) 



For equation (A5), integration yields 
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(48) 



and for equation (46) 



(49) 


With these intermediate results, the shape parameter 6*/G can 
be given in a form that now only depends on the corresponding 
temperatures, and T^. 


Q 


2 

c g_ 




T T- 0.5 

" 2 ♦ 

c eg 


■ t 

D.T 


♦ 1 


- 1 


(50) 


Figure 14 (page 39) shows a plot of this parameter for increas- 
ing Mach numbers (or e, respectively), with y as the parameter. 


For a non-adiabatic wall, the wall temperature will be lower 
than the stagnation temperature, and the gas density near the 
wall will be higher than in the adiabatic cases. Therefor 6 * 
will be smaller and also 6*/0, because 
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‘ 2 5 10 20 _ 50 100 150 250 

Surface ratio e 

Figure 14 Increase in boundary layer shape parameter 6*/© as a function 
of e with Y parameter 



The gradient 6*/0 is not much affected, however, precisely for 
the case of large surface ratios, in which the stationary wall 
temperature is considerably below either the combustion chamber 
or the stagnation temperature. 

2. 2. 5. 3 Nozzle shape correction 

Due to the formation of the boundary layer, the inviscid geo- 
metrical cross-sectional area of the nozzle is reduced by the 
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impulse loss thickness. To determine the resulting thrust loss 
it is necessary to determine the fully viscous velocity distri- /3A 
bution, which can be derived assuming that the mass throughput 
through the narrowest cross-section also flows through the 
nozzle end cross-section (contiguity equation) 


ID • m « ■ n 
V n/v 


(52) 


Including the usual boundary layer quantities, we obtain: for 
the mass throughput 


#0 /f* 

g.v ^g.v 


♦ (r^ - 


g.v g.v 







(53) 


and for the thrust, 


F » TI . U* • Pg (r^ - 6 )^ ♦ 2tI (r^ * 7 ) * Pg * “g * * 


J p U* dy 


(54) 


with the abbreviations 
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g.v 


•^g.v 

p/p, 


u/u 


g.v 

g.v 


= p 

■ p 

« U 


We can now rewrite equations (53) and (5A), as 


itu 


p„ (r - (r - |) P,u,« - r> ' Vv ■ %/v 

40 


V 


HU. 


(r, - 6)' 


♦ 2ii 


(r, - |) 


0 u 

V 


fi «* 

-r 


(56) 


The viscous thrust {includes: boundary layer behavior due to 
gas viscosity and wall friction) will be, expressed in abbrevia- 
tions: 


'n/v 


Vv • 



(67) 


Using the customary definition for specific impulse: = 

F/m we obtain from equation (57), taking into consideration 
equations (55) and (56), the following loss of impulse 


AI 


sp.bl * ^0 “ “n/v ' u 


7 V 7 “ 2 

n/v n/v 


P 

V » 

“n/v'%/v 




(58) 
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Al *0 
sp ^0 


Vv 


2 

“n/v*®n/v 



(59) 


The denominator of the factor in the square bracket (Kg) can 
also be transformed, with the aid of equation (55), in a form 
that depends only on the boundary layer parameters 


^n/¥ 


^n/v 


’>*Uy’Py(r-6)^ 


2h (r - ^)p 


2'W 


(1-4^) 


V r* 


(60) 


or respectively, 
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“n/v * “v 


Equation (59) for the loss of impulse now becomes 


sp ’o n/v "v IT 

n 


The determination of the velocity ratio 


^ 1 
Vv “* 


occurs by means of the relationship already presented in equa- 
tion (30), and with the known definition of the Mach number 


m 2 “g 

'•s'VTsV 


To begin with, we shall divide our considerations into two 


cases 


■ Vv • ' 




Case 2 


Case 1 ) : 


. 1 ♦ Y'l 


1 . Y-1 Vv 


With 


1^ - ^ Y-1 ,-0.5 

\ ‘mn^' 


theref ore (66) 


Vv " n - 

c c 
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The abbreviation Kj is a constant for local nozzle points. 

c 

A similar treatment follows for Case 2): 


u 


V 




0,5 


( 68 ) 
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Now the velocity ratio in equation (63) can be expressed in a 
form that depends only on T 

< 1 ^/ ■ - T„/,) (69) 

With equation (61) and the general gas law, as well as a rewritten 
equation (69) 


(70) 



therefore 


_I!:L . 1- JiL (71) 

“n/v ^n/v 

Since the pressure can be determined from a temperature function 


v-1 

p y = T 


y 

p » T 


(72) 


we conclude 





( 


^n/VjY=r 
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(73) 


But since the viscous temperature is not known explicitly, 
is replaced by equation (70) and equation (73) now becomes 



or respectively. 


it 


( 74 ) 




<• 4 , • “ 5 > • “ 4 , 


“S/,’ 


Y*1 

Vv 


• A 


(75) 


The right hand side of this equation is already known, through 
the quantities that can be determined by the potential theory 
and hence the fully viscous velocity u^ sought must be 
determined by iteration from the final equation, which can not 
be solved explicitly, 





(76) 


By means of the very important quantity u^ here determined, it 
is now possible to incorporate the influence of boundary layer 
development on the core flow within the nozzle, directly into 
the computer program. 


2.2,b.k Correction of the compressive thrust term 

In analogy to the influence by the viscous velocity development, 
the compressive thrust portion must be corrected to include the 
boundary layer conditions. 


The entire loss of thrust due to boundary layer effects is 
expressed by 

^bl * * ^e,2-d ^Pe.n/v * ^e.v^ ^77) 

Taking into account the results of the u^ = f(u^^^) determina- 
tion, the following relationship exists between the pressure 
equation and the velocity distribution 


u 


\ - ( ^e.n/v .y 
Pc 


1 - 


P ^ 


( 78 ) 
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Rearranging and solving by viscous and inviscid terms yields 


1 


(Vv,2 

u 

V 


P ; 1=1 

^ c.n/v ^ Y 




P ^ 


•■s. 


( 79 ) 


The right-hand side includes only magnitudes already determined /^1 
and hence the fully viscous end pressure now depends only 

on Pg, y and the already determined value Kp. 


^e,v 


Pc S 


( 80 ) 


2. 2. 5. 5 Connection between nozzle end values and boundary layer 
losses 


Based on the preceding considerations and the Fortran IV computer 
programs built on them, all required fully viscous nozzle end 
values can be numerically obtained with sufficient precision. On 
this basis, the total integrated boundary layer loss for super- 
sonic flow can be obtained 


81 


sp 


• 8o 


u 

“n/v,e ' u 


2 

v.e 

n/v,e 






• (1 - (1 - -jr - r-) ♦ Ag - p^ 



( 81 ) 


or, related to a 1 -dimensional loss free theoretical impulse 




^sp,n/v ~~^^sp.bl 
^sp,n/v 


( 82 ) 


A review of the numerical treatment of the functions and differ- /A2 
ential equations in the preceding portion is contained in the 
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form of a flow diagram, in the Appendix. The integration of the 
impulse differential equation (see section 2. 2. 5. 2) is performed 
applying an explicit differences procedure, based on the Lax- 
Wendroff method [35 ,8A ,9^ , 107 i 115] . 

2. 2. 5. 6 Polynomial for the velocity profile 

In order to prove that the simplified linear theory of the 
velocity profile described above does not represent an impermis- 
sible imprecision, a polynomial profile will be used in addition, 
derived from laminar flow (Pohlhausen) [136,1653, 

^ • (2n - • n^) ♦ J (n - , 3n’ - n‘) (M) 

A * ^ ^ 0<y<6 

The results thus obtained show no decisive departures from 

the results obtained with the linear profiles and hence all sub- 
sequent cases will be investigated using it. 

2. 2. 5. 7 Roughness equation 

To include the effect of the roughness of the nozzle’s inner 
wall on the formation of the boundary layer, at each computation 
stage in the core flow, Re^^^will be compared to the values 
known for pipe flow. As is known, at Reynolds numbers >2,300, /A3 

Colebrook-White ’ s equation [161] is valid for the transition 
range of the pipe friction number 

1 _ o T ^ „ 2, SI f, Rdu 

In contrast, by calculation of the laminarization parameter it 
is now known that laminar flow predominates. Hence, the calculated 
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Reynolds number {> 2,300) must be recalculated - with the aid 
of the pipe friction coefficients - to the Reynolds numbers of 
a laminar pipe flow which is possible with 

Re, . l«) 


The effect of roughness is seen by correction of the viscosity 
parameter 



Re 




n/v 


(87) 


In the calculation processes, wall friction is taken into con- 
sideration by means of 

(k = 0.5 in the smooth case). 


In this model the factor k = f(rough) is connected with a 
Reynolds number that takes into account the displacement thick- 
ness 6» and indirectly - through the corresponding surface ratio 
e - also includes the geometry. If we set k = f(rough/R), then 
the roughness is correlated directly with the geometry (Figures 
15 and 16), as generally customary in pipe friction theory. 


Other possibilities for a roughness equation are listed below. 

A roughness expression in terms of the wall shearing stress coef 
ficient c^ seems most appropriate, since it is affected by the 
roughness of the surface layer. Thus, all possibilities that can 
expand by means of an additive or multiplicative term, are 
open . 


An additive term 


c . 0.5 



A7 


/kk 


(89) 




Figure 15 Effect of rouguness calculated using REKOR. Nozzle 5 

offers the possibility of introducing = f(rough) and to form 
the Reynolds number Re^ with the core flow velocity, or the velo- 
city at the height of maximum roughness; Reynolds numbers formed 

on the basis of mean velocity values are also conceivable. A 

0 5 

multiplicative term c- = — k- makes it possible to define a 

kind of loss efficiency of the wall roughness, in which the 
height of the roughness is related to length, to form k^. Suit- 
able lengths are 6*, r^ or a combination of both. The literature 
does not offer any suitable roughness expression that would be 
universally applicable to small nozzles with Re > 2300 (laminar). 
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(M6) 



RAU (rough) in m 


Figure 16 Effect of roughness calculated using pipe friction expression 
(REKOR) 


Hence, the best suited must be found by measurements [15], 

Anticipating the results of the tests here described, we find 
that roughness has no substantial effect on » which confirms /A5 
the approximate expression used earlier. The explanation for 
this could be that the relatively large boundary layer thick- 
nesses in the case of small diameters extend far beyond the 
maximum roughness heights, such that the latter are unable to 
develop any large loss effects, i.e., the roughness effects come 
to rest within the viscous boundary layer. (Microscopic measure- 
ments of the manufactured propulsion units show an absolute 
roughness of < 1/1000 mm). 

2.2.6 Adaptation losses due to environmental pressure effects /k'j 

For a given pressure distribution within the propulsion unit, 
the thrust can be determined by means of 
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Integrating over the entire inner surface of the engine » we 
obtain, in a scalar decription for 1 -dimensional flow, 

F • roc ♦ (Pg - p.) Ag (91) 

The first term represents the so-called impulse thrust, while 
the second term, formed by the pressure difference between the 
nozzle end and the environmental pressure, determines the 
compressive thrust portion. 

Since the expansion ratio of a nozzle is not infinite at any 
time, Pg always has a finite value greater than 0. Since in ad- 
dition in the actual test case p^^^ is greater than 0, the propul- 
sion unit's measured thrust must be corrected to attain the 
vacuum thrust necessary for performance specific comparisons. An 
additional precondition to this is that during experimental 
operation fully formed nozzle flow is maintained, without 
separation phenomena and compression shock within the nozzle, 
such as those occurring in the case of superexpanded propulsion 
units p << p„. 

(For large nozzles, the Sumraerfeld criterion can be applied to 
good approximation: p^ > 0.4 p^„; in the present case this ratio 
is shifted downwards - see also section 4.2). Since in reality 
the pressure distribution across the end cross-section A^ is 
not constant, a 2-dimensional pressure distribution is used to 
correct the measured F values, based on the calculation method 
and the theory described in connection with section 2.2.5. 

2.2.7 Heat losses 

In contrast to the large chemical propulsion units, which show 
a substantial decrease in wall temperature along the axial 
direction, here a homogeneous temperature distribution is 
assumed. The following boundary conditions must be mentioned, 
in this context: 
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a) Large propulsion unit heat capacity 

b) Large heat capacity of suspension and supply lines 

c) Short test duration 

d) Good insulation of the entire assembly in the vacuum 
chamber 

e) Small temperature differences between thrust chamber 
and environment 

f) Radiation loss remains small, due to 

As a consequence of the last three points mentioned, heat losses 
can be neglected, in the comparative calculations for propulsion 
unit performance. Performance measurements already carried out, 
in conjunction with high-resolution temperature sensors for the 
control of adiabaticity , show this to be permissible. 

2.2.8 Results of performance prediction calculations 

For all the propulsion units examined in the course of the experi- 
ments, the most important thermodynami and fluid dynamics para- 
meters were determined using. 'the EGVP procedure introduced. For 
better visibility (and also due to efficiency considerations), 
these parameters were plotted as a function of e, using a specially 
develeoped boundary layer program (or alternatively, as a function 
of the nozzle axis coordinate). 

Figure 17 (page 52) shows the pressure, temperature, viscosity (/50) 
and density for nozzle No. 2, as an example. In addition, both in 
this and the two subsequent diagrams the nozzle shape was also 
drawn in. 

/ 

Figure 18 (page 53) represents the changes in the core flow (/5l) 

velocity, Mach number and larainarization parameter, as a function 
of the no-zzle length z [A5,^6]. 

The courses of the boundary layer parameters shown in Figure 19 
(page 5A ) show that the boundary layer thickness in the throat 
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Figure 17 Course of the pressure, temperature, viscosity and density 
in nozzle No. 2 


has been reduced to practically zero, in the convergent portion 
of the nozzle, by the predominating strong acceleration of the 
gases in the previous range (see also section 2. 2. 5. 2. 3). 

With the two velocity profile expressions introduced in section 
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Figure 18 Course of fne core flow velocity, llach nujsber and laminar izat ion 
parameter K along the nozzle axic (Nozzle No. 2) 


2.2.5 and with as parameter, Figure 20 (page 55) shows typical 
boundary layer losses, using nozzle No. 2 as an example. It can 
be seen that the linear and polynomial velocity expressions have 
no qualitative effect on the results. As p^ increase, there is 
a considerable decrease in losses. The quantitative differ- 

ence (a factor of approximately 2) is based only on the expression 
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of A and can be readily adapted (see also section 2. 2. 5. 6). 

Figure 21 (page 56) provides a comparative representation of 
the losses various media. Qualitatively, the boundary 

layer losses as a function of c can be interpreted as reflecting 


/53 


CMM) 



Figure 20 Boundary layer losses without roughness correction 

the effect of the different viscosities (see also section 
2,2.10). The only exception to the general trend (increase of 
^^bl increasing e) is that of the curve for (NpH, ) 

^dis = 1.0 

Due to the strong cooling of the originally fully dissociated 
ammonia, a recombination to NH^ is initiated, with the simul- 
taneous release of energy, allowing the gas better expansion 
performance. Consequently the bosses decrease, percentually . 

In contrast to large propulsion units (F-1, J-2, SSME), an addi- 
tional design optimization is necessary with small propulsion 
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Figure 21 Surface ratio as a function of boundary layenloss for various 
types of gas., = 15* 

units, because of the greater influence of the boundary layer 
losses: the optimization of the surface ratio e with respect 

to and (not to be confused with the e-optimization of 

carrier rocket propulsion units with respect to p^ r f(H) and 

H)). In general the losses increase with 

increasing a„ • Since the specific impulse I steadily 

increases with increasing e, an optimization calculation has to 
be performed to obtain the optimum e. This must be followed by 
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Figure 22 Boundary layer loss effect on the determination of the 
optimum surface ratio (N 2 H 4 , rj- = 0.3 nnn, = 15 *) 


a second optimization, to take the influence of into account: 
while at an assumed constant e, for an increase in a^, becomes 

larger (conical nozzle), at the same time the absolute cone length 
decreases and hence the development of the boundary layer is fore- 
shortened. The mathematical formulation is 


I 


sp, ODE, vac 


1 

^ “ TOT j W 


I 


sp,max,vac 


( 92 ) 


or 


j^^p.ODE.vac * * TOT )]“ ° 

Figure 23 {page 58) shows the result of the study for the 

nozzle, as an example (56 t 63). 
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Divergence angle a^, in degrees 


Figure 23 Optimum divergence angle taking the boundary layer and multi- 
dimensional losses into consideration, with £ as parameter (N 2 H^) 

Since the optimized I , , curve shows a very flat maximum at 

S ^ ^ O 6 

e = 80 approximately, and in addition the effect of various 
errors is restrictive (for instance, the true degree of dissocia- 
tion and the recombination behavior as a function of the expansion 
and the residence time in the TW [Triebwerk = propulsion unit], 
a fairly large bandwidth must be assumed here. 

Figures 2k and 25 (pages 59 and 60) show the optimization /56 
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results for nitroher. Is addition to the line connecting the 
corresponding according to the explicit theory, 

Figure 24, above, also shows the line obtained connecting the 
points resulting according to the simple approximation equation 
(see section 2.2.10). The good agreement is apparent. The 
computation time decreases by three orders of magnitude, in the 
process. Figure 25 (page 60) shows the resulting values, 

as a function of r^, for a constant value of Due to the low 
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viscosity of N 2 iri comparison to, for instance) H 2 or N 2 H^~ 
decomposition products, the optimized surfscQ ratio becomes 
displaced towards larger values. 


2 . 2.9 Comparison of the results presented in section 2.2.8 with 
those from other procedures (Cline, BLIMP [Boundary Layer 
Integral Matrix Procedure]) 

The flow mechanics quantities calculated by the EGVP program 
- pressure Pg, temperature Tg, velocity Ug and Mach number M - 
are compared to those from the program VNAP (Cline [22,23,153]), 
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a program for 2-dimensional, compressible, time-dependent, fully 
viscous nozzle flow. Figures 26 to 29 (below and pages 62-64) 



Figure 26 Comparison of programs EGVP and VNAP with respect to T for 

1^2 nozzles 

show the results. In the calculations with VNAP, the values of 
the core flow present discontinuities just before the nozzle 
throat; at high surface ratios they coincide with the values 
calculated by EGVP, which also shows good agreement in the 
throat region [ 50 ]. 

We have no explanation for the oscillations before the throat, in 
the values calculated by VNAP. The departures behind the throat 
could conceivably stem from the fact that VNAP uses only a throat 
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Figure 27 Comparison of the programs EGVP and VNAP with respect to the 

course followed by the pressure 


transition radius and hence is unable to describe the optimum 
nozzle geometry with sufficient precision [70,102,105,106]. 


Since both the EGVP and the VNAP programs determine fully viscous 
flow mechanics data, and since derived from them, we 
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Figure 28 Comparison of the programs EGVP and VNAP with respect to the 
core flow velocity u 

g,v 

conclude from the agreement of u^, M, and found above, that 
the boundary layer development is comprehended with a similar 
degree of agreement. 

While work with the Cline procedure seemed relatively simple, 
due to correct program listings and a very precise description 
and good collaboration with the author, this can unfortunately 
not be extended to the second main comparison program, BLIMP 
[ 16 , 113 ]. 

Through at first fortunate circumstance, the Boundary Layer 
Integral Matric Procedure (BLIMP) reached the Chair for Space 
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Figure 29 Program comparison EGVP-VNAP with respect to the 

Mach number 

Technology. This program had been developed during the Seventies 
in the US, involving large expenditures in materials and leading /59 
specialists, in the context of JANNAF, for application to many 
gas dynamics problems: from the flow over a flat surface under 
subsonic conditions, to 3-dimensional, non-stationary supersonic 
flow in rocket propulsion units. Since the documentation is 

6A 
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Figure 30 Course of boundary layer paraioeters obtained using BLIMP 
p^=10 bar; T^=2300 K; r^=0.2 mm; a^=15®; 

insufficient and additional, machine-specific problems arose at 
the Leibnitz Computer Center of the Bavarian Academy of Sciences, 
we have to date been unsuccessful in our attempts to use the 
BLIMP program at the Technical University Munich. The example 
shown in Figure 30, above, was obtained directly in the US. It 
shows the boundary layer quantities 6, 0 and 6^ (energy thick- 
ness) for a combustion chamber pressure of 70 bar, T^ = 2,300 K, 
r^ = 0.2 mm, = 15° and (decomposed and heated). With the 

help of the energy thickness 6^^ (see for instance Truckenbrodt 
[167]), it is possible to provide more comprehensive loss 
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statements than with 6*, otherwise {EGVP) used here. The 
mathematical formulation follows a similar principle and causes 
no additional problems. The comparison of BLIMP and EHVP data 
in regard to propulsion units shows sufficient agreement, 

and hence additional justification for the simplified EGVP 
procedure and the system developed in section 2.2.10 is not 
necessary. It must still be remarked that a BLIMP calculation 
requires an initial data set provided with the aid of ODE, ODK 
(Gordon McBride) and TDK. After the first run of the BLIMP sys- 
tem the nozzle shape is corrected and the procedure is repeated. 
Customarily satisfactory convergence is established after two 
passes. Depending on the geometry and the reaction chemistry, 
the computer time required varies from 30 min to 2h CPU time. 

For comparison purposes, EGVP requires from 5 to 30 seconds. 

2.2.10 Simplified performance prediction model 

While the computer program for predicting the performance of 
small rocket propulsion units, described in the preceding section 
is considerably more modest, in regard to computer time and 
expense (by a factor of 100, compared to BLIMP), it still does 
not satisfy the criteria of simple engineering "rule of thumb" 
equations. This shall be attempted in the example of boundary 
layer losses, below. 


2.2.10.1 Derivation of the approximation equation 

The starting point is also the impulse equation, slightly 
modified for 1-diraensional flow, 


do 0 tan a 3 ♦ 2 - M' 

Si “ Z"cos~ * ^ J.s* 

■ 


(94) 


using the Reynolds number related to the impulse loss thickness 


Re « 


PcO 

c* c n 


( 95 ) 


and with the viscosity function 
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. k , 


( 96 ) 


(The e-dependence occurs because of the application of the ODE 
theory to T^/T). 

The shearing stress coefficient is determined as 

Uk 


^S.c ^ 

c. ail 


^n,2 


p,0 


'f.2 


(97) 
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The impulse equation can now be rewritten 

d0 _ 11 ^ 




15 '^G.1 


m 2 

The surface ratio e is replaced by (^) . Using 0 = 0 

(start of integration at the narrowest cross-section) , we obtain 


(98) 


(99) 


filial k r 

•'o' t 


or 


( 100 ) 


We can now determine as 


A- "5.cV'* 


‘‘n.bl.l (^*«Ve*6‘/2r^O) (101) 


Introducing the throat Reynolds number 
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the expression for the boundary layer losses can be strongly 
compressed : 
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2.2.10.2 Loss estimates for various working fluids 

With increasing Reynolds numbers in the throat region - which 
can mean higher chamber pressure, larger propulsion unit dimen- 
sions and higher thrust - the boundary layer losses decrease, in 
terms of percentage. In contrast, an increase in the surface 
ratio also increases the losses. It is important to establish 
that the trend of equation (103) remains unchanged, even using 
other expressions for c^. This has made it possible to formulate 
available experimental data and the theoretical results of 
extremely comprehensive computer programs with the help of a very 
simple equation. Hence the design of a future propulsion unit 
can be determined very efficiently with respect to the optimum 
surface ratio e. The precision of the result of equation (103) 
is surely sufficient, if we take into consideration how difficult 

it is to determine the other parameters, such as c* , p , T , e, 

c c 

etc. To determine the two constants in equation (103), a series 
of available experimental data and voluminous calculations were 
performed with the corresponding procedures, for different media. 

The results of the determination of correlation coefficients for 
various working fluids are shown in Figure 31 (page 69). The 
Freon- 1 A curve lies above those for simple gases, even though 
the viscosity of for instance, is higher than that of CF^. 

Apparently, the explanation for this departure must be sought in 
that during the expansion process (and the concomitant cooling 
of the gases) CF^ recondenses to small droplets (boiling point 
at 1 bar approximately -130°C), leading to higher loss measure- 
ments. 
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Figure 31 Course of the boundary layer loss coefficients for various 
types of gases, from available experimental data 

In an ideal case, the curves shown in Figure 31 should come to- /67 
gether at one point, for each type of gas. If this does not occur 
it is because the diagram was built using many sources, whose 
individual results show considerable departures [60] . Because of 
the hyperbolic nature of the curves and the mathematical quality 
of equation (103) it is possible to insert pairs of values 

(by prescribing e and , without the boundary layer loss 

diverging substantially. 
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3 EXPERIMENTS WITH MICROPROPULSION UNITS TO BETTER CORRELATE /68 
THE RESULTS OBTAINED IN SECTION 2 

3.1 Requirements for the test and measurements facilities 

A substantial experimental program is carried out to provide 
an underpinning for the results attained in the preceding section, 
to measure specific performance quantities in miniature propulsion 
units. As can be anticipated from the theoretical results, severe 
requirements must be placed on test design and measurement preci- 
sion, in order to be able to identify the individual performance 
losses, for instance due to boundary layer effects. 

The development of detailed test facilities must start, in each 
case, from the following boundary conditions and physical 
quantities : 

TABLE 3 Precision requirements for the measurement installation 

Thrust range: 10-2000 mN ± 1% 

Mass throughput: 10-2000 mg/sec ± 1% 

Combustion chamber pressure: 1 - 20 bar ± 0.3% 

Combustion chamber temperature: 290-300 K ± 0.5 K 

Nozzle throat diameter: 0.25-0.8 mm ± 0.003 mm 

Surface ratio: 1 - 200 ± 0.5% 

Gas type: repurified (99.99% N 2 ) 

Measurement durations: 1 - 100 sec ± 0.1% 

Safety 

Reproducibility 

Manufacture of individual components in the Chair's own shops 
Financial context < 10,000 DM (outside financing) 

3.2 Test facility description 

Without the Chair for Thermodynamics (Prof. Dr. E. R. F. Winter) 
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Figure 32 General view of the altitude test stand for cold gas propulsion 

units 

at TUM [Technical University at Munich] generously providing a 
laboratory, due the extreme space scarcity at the Chair for 
Space Technology this experimental study in the area of rocket 
propulsion unit performance measurements would not have been 
possible [55,163]. A massive frame construction was installed in 
the already available (as a loan from the DFVLR [Deutsche Forsch- 
ungs- und Versuchsanstalt fuer Luft und Raumfahrt r German Research 
and Test Facility for Air and Space Flight]) altitude simulation 
chamber. This frame holds the propulsion unit's suspension, F- 
calibration device, gas supply, p^ and p^^^ measurement pipes, 
temperature sensors and components of the optical thrust measure- 
ment arrangement. The Figures 32 (above) and 33 (page 72) show 
an overview. The N 2 supply storage (including the m measurement 
system) is located outside the vacuum chamber, as are the main 
components of the thrust measurement system, various electrical 
supply devices and measurement value acquisition and processing 
[ 17,20,40,52,7^,75,78,98, 101 , 155, 168, 169, 178] . 
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EXY to Figure 33 (page 72) 1 Manometer (quality rating 0.6) 2 Exhaust 

gases 3 Timer 4 Needle valve 5 Electromagnetic stop-valve 6 2-stage 
blocking piston valve piunp 7 Printer 8 Main valve (gross adjustment) 9 
Electromagnetic bypass-valve 10 lO-channel input 1 1 Multiplexer 12 A/D 
converted 13 Memory 14 Control 15 Printer 16 PD regulator 17 Kombitron 


3.2.1 Vacuum installation /71 

To compensate for the upthrust of the thrust equipment during 
calibration of the thrust meaurement system, the environmental 
pressure in the vacuum chamber (useful volume approximately 150 
liters) must be adjusted to the underpressure present during the 
calibration phase of the ’’nozzle test run” to be performed sub- 
sequently. To this end, the suction of the rotary vane pump (LH- 
S60) can be adjusted by means of a ’’bypass” needle valve. A 
specially developed analog PD regulator (see Figure 34, below) 


Nominal value 



controls the ”opening” and ’’closing”, respectively, of the bypass 
adm.ission of external air into the rotating vane pump’s suction 
line. By using LH-Thermovac and LH-Penningvac pressure measure- 
ment pipes in connection with the LH Combitron supply and 
indicator system, the actual value of p^ is supplied to the 
control circuit (Mr. Matthern of DFVLR Oberpfaffen was kind enough 
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to perform the calibration of the sensor system for us). 


The Pjp adjustment makes it possible to achieve a reproducibility 
of better than 1% in the pressure range of 10“^ - 10 Torr here 
under consideration , which is sufficient with respect to thrust 
measurement errors due to upthrust effects (see also sections 
3. 2. 2. 2 and A.3) . 


3.2.2 Thrust measurement system 
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In general the thrust of a rocket propulsion system can be 
described as the product of the effective discharge velocity and 
the mass throughput: 


in c. 


eff 


( 104 ) 


While in most test cases the propulsion unit’s mass throughput 
can be determined with reasonable efforts, a direct, precise 
measurement of is not possible. This is due to the follow- 

ing reasons: 


- 3-dimensional distribution of the discharge velocity 

- 3-dimensional pressure distribution in the end cross- 
section 

- Condensation shock, flow separation 

- Assignment of local density and velocity distribution 

- Sensor problems 

- Reduced size of the nozzle geometry (in the case to be 
examined here) 

Hence only an integral measurement of F can be used here. Because 
of the small forces involved, the methods used for large rocket 
propulsion units - using DMS sensor to measure the deformation 
of the thrust stand - can not be applied. While inductive 
displacement pick-ups could measure some of the "large” nozzles 
(2N) used here, their resolution is limited for the lower thrust 



range (20 mN), and they can hence not be used. It was decided 
not to use torsion balances - as they are used for small 
single-component fuel engines (ERNO) - since they aloow the 
determination of only the thrust time integral, from the total 
torsion angle, rather than the temporal course of the thrust. 

Because of the drift problems in the charge amplifiers that would 
necessary for the nearly stationary operation case, piezoelectric 
sensors can not be used. In addition, to reduce the mechanical in- 
fluence of supply lines and calibration devices on the thrust /73 

measurement method, contact-free methods should be included in 
a narrower selection. 

- Due to the required resolution, it is not possible to 
use radar methods 

- Optical interference methods are possible, but very 
demanding in their implementation and difficult in terms 
of measurements evaluation (linearity errors, counting 
problems ) 

- Optical "edge" and "slit" system were evaluated in pre- 
liminary testing but eliminated because of their strong 
dependence on the light source's intensity value. 

- The situation is analogous for systems that translate 
the smallest deflection of thrust transverse beams into 
intensity changes, by means of the displacement of a 
polarization filter 

- The F-measurement principle finally chosen is described 
in detail in Figure 35 and 36 (pages 76, and 77): 

A simply polarized laser beam strikes a mirror installed on the 
propulsion unit's flange, parallel to the thrust direction (first 
surface mirror) . The thrust frame consists of a cantilever beam 
that simultaneously represents the gas supply pipe. At its free 
end is the installation flange with the propulsion unit, which 
generates a thrust acting vertically downwards. This causes the 
pipe to be deformed along its deflection curve bending line. 
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such that the mirror alters the angle to the incident and the 
reflected beam. Outside of the vacuum chamber the reflected beam 
is sent to a homogeneous linear array (LSC-5D, United Technol- 
ogies). This beam position measurement probe - which operates on 
the Schottky principle - permits a resolution of better than 
0.1 pm, using the appropriately adapted electronics { TS-Elektronik 
Munich, FRG ) . 

In spite of the .required resistance of the bending pipe against 775 
internal pressure (to a maximum of 25 bar), the resolutions allows 
dimensioning that pipe in such a manner that in conjunction with 
the corresponding optical "lever arm" thrust forces in the range 
of 20 mN to 2 N can still be measured with sufficient precision. 
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Figure 36 Mechanical test construction within the vacuum chamber (shown iso- 
metrically and partially diassembled for better understanding) 


Interpretation involves purely elastic resistance expressions, 
which after numerous calibration measurements seems to be permis- 
sible (se also section 9.1). 


Since both the LSC array sensors and the necessary evaluation 
electronics still require empirical values from use in day-to-day 
laboratory practice - in contrast to the segmented position sensors 
already used for longer periods of time - it may be assumed that 
the reproducible resolution of the position measurement procedure 
can still be improved upon- It does not seem necessary to point 
out in detail that the method here presented is also suitable for 








larger propulsion units and other problems, in which the record- 
ing of very small position changes is necessary. Although in 
principle the optical thrust measurement procedure is of very 
simple design, during implementation the edge effects causes 
problems whose solution is quite laborious .Diffraction and refrac- 
tion influences of the laser beam, caused by optical windows and 
beam refraction on gas molecules, can be corrected, or even avoided 
by the use of polarization filters, quartz windows and the required 
air density control in the vacuum chamber. 

The main problem occurring in the intended experiments is keeping 
the laser beam intensity constant. 
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Figure 37 Laser-optical thrustmeasurement system (beam path 
made visible using smoke. Exposure tine: 4 minites) 
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3.2.2. 1 Laser intensity control 

Long and short-period variations in the laser beam output 
intensity are due to line voltage variations (long-term and 
short-term) , temperature influences and the plasma physics 
processes inside the laser tube. In the case of the He-Ne laser 
system used here, the amplitude can be up to 2% of the nominal 
initial output. In addition, there is the constant degradation 
of the beam output due to diffusion of foreign atoms in the gas 
laser's glass tube. The array diode deflection measurement 
procedure is capable of linearizing and compensating for the 
changes mentioned above, in a first approximation (including 
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the subsequent amplification and evaluation electronics). The 
remaining higher order interference terms, however, prove to be 
intolerable in their effect on the resolution and reproducibility 
of the thrust measurement procedure. 

Hence it is necessary to develop a procedure for an active 
stabilization of the laser intensity. To this end, we first 
replaced the laser system CWR-LS-05 (with unregulated power 
supply), by means of the model Laser Optronics LO-101 P, whose 
power supply reduced the initially necessary ignition voltage of 
10 kV, to a continued operating voltage of 1500 V, after the 
start of operation. In addition, it smoothes out intermediate- 
term tension variations of the supply line voltage. 

We decided against the purchase of a very expensive (nearly 
DM 30,000) power supply, especially tailored for laser problems, 
since it can eliminate only certain line influences in the short 
and long-period fluctuation range, but does not eliminate the 
thermal variations of the laser intensity, or those due to plasma 
physics, ultimately. Among them are, for instance, the effects 
caused by "mistuning” of the laser mirror system at the two tube 
ends, caused by thermal shifts. 

The intensity stabilization system finally used was shown in /78 
Figure 35, in schematic form. Via a prismatic beam splitter, 50% 
of the beam output is deflected towards a photodiode operating 
in the He-Ne band, as an ’'actual value" indicator (in order to 
compensate for the influence of the different gas densities on 
the diffraction of the laser beam, the beam splitter is arranged 
within the altitude simulation chamber). In a specially developed 
PID regulator (circuit diagram below, page 81), the instantaneous 
actual value of the laser intensity is compared to an adjustable 
nominal value [AO]. When the appropriate adjustment difference 
occurs (= 50 yV) it is amplified sufficiently so that a micro- 
metric control can be set in motion via a down-stream, direct 
current minimotor. On this mechanically prestressed control (to 
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Figure 39 PID Laser intensity regulation 


avoid slippage during pole reversal), there is located a 
continuous neutral wedge. This optical component - which acts 
as a linear filter - is placed in the beam path, between the 
laser exit and the prism mentioned above. It thereby becomes 
possible to correct for the smallest intensity variations and 
the continuous output decrease (life-tiem dependent). The integral 
portion of the regulator limits the long-term weakening to less 
than 3%. The differential portion of the control electronics 
makes it possible for the regulating circuit to react already at 
the beginning of a steep interference flank, before maximum ampli- 
tude is attained (leading effect). Large demands are placed on 
the proportional part, because of the necessary amplification 

factor of approximately 10^, since circuits of these character- 

3 

istics, with amplification factors >10 tend to be dynamically 
unstable. Hence, a multi-stage amplification chain should be 
adopted. It is absolutely necessary, in addition, to use 






"offset"-balanced operational amplifiers of a model with high 
temperature stability specifications. In spite of this very 
laborious detail work, in its technological implementation, we 
were at first not sully successful in compensating for the very 
strong short-term line fluctuations, determined by specific 
local conditions (caused by the discontinuous operation of 
machines of high power consumption). It thus became necessary 
to operate an independent supply line (see Figure 52). It consists /79 
of a 2A V, 15 A battery charger, with a AOO Ah fork-lift 
accumulator and a 24 V DC/220 V, 50 Hz direct current/alternate 
current inverter (200 W output power). A test/loading cycle time 
of approximately 1:1 is thereby attained. The nominal block to 
block time is approximately 20 h: after at the most this period 
of time the laser power supply must be switched back to the 
public power network, to recharge the accumulator. Since the 
entire laser-optical thrust measurement system is kept under a 
light-tight cover to avoid outside light influences, active cool- 
ing must be provided to prevent heat accumulation near the laser, 
due to the laser’s own heat losses. Open water circulation circuits 
are adequate for this purpose. 



Figure 40 Portion of the thrust measurement system outside of the 

vacuum chamber 
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Figure 41 Neutral wedge control (active control member) for laser 

- intensity stabilization r-r 

3. 2. 2. 2 Calibration device for thrust measurements : /81 

In order to calibrate the thrust measurement system, first the 
environmental pressure occurring inside the vacuum chamber •: 
is established, in a preliminary tests, as it is under full ,1: 
pump suction power and nozzle flow (leak rate) at fixed p . Next- • 
the N„ supply pipeline is interrupted and the p regulating ; 

circuit is activated (see also 3.2.1). Thereby the environmental: 
pressure is fixed at the value determined above and the upthrust-^ 
influences due to residual gases in the chamber may be considere:d 
compensated. Now the calibration basket (see Figure 43, page 84,):. 
is introduced through the divergent nozzle portion. The basket -s: 
positioning is so designed (double cone, point contact) that 
calibration masses placed in the calibration basket can exert 
only vertical components in the thrust direction and position of, 
the nozzle axis. By means of a revolver cassette system, arbitrary 
quantities of beads of precisely know mass can be admitted to the 
basket. These beads simulate the "calibration thust" for the 
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Figure 42 Qualitative course of the thrust calibration 
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Figure 43 Propulsion unit suspension with flanged combustion chamber 
pressure determination and calibration device of the thrust 
measurement system 


thrust measurement system. The response of the measurement system 
is shown in Figure 42, above, for the individual beads falling 
into the basket (Mass: typically 3 g). 
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From the step function obtained during the calibration phase /82 
it is now possible to determine the gradient of the function 
F = f(U Once the calibration measurement is com- 

pleted, the calibration basket including its position cone is 
removed from the propulsion unit by means of automatic controls, 
so that that the precision performance measurements can be 
performed immediately following. Unavoidable zero-drift in the 
instruments is thereby prevented from exerting undue influence 
on the precision of the measurements. For additional control, 
the calibration was repeated after each measurement series. 

3.2.3 Mass throughput determination /83 

The mass throughput m is the most important limiting quantity 
- next to the thrust - in the experimental determination of a 
rocket propulsion unit’s specific performance Igp* 

A large number of physical principles - which in turn can be 
expressed technologically in an enormous multiplicity of forms - 
could be used in its measurement. In the context of this study, 
the most significant possibilities were examined either by 
preliminary theoretical considerations, by borrowing existing 
systems and by extensive preliminary testing. Therefore, the 
listing following should be consiered only as an overall review. 

3.2.3. 1 Comparison of the various alternatives [28,31] 

a) Throughput measurements using turbine rotor systems: 

A turbine rotor is fitted into a feed pipe, with its segmented 
axle made alternatingly of conducting and non-conducting 
material. The fluid's motion causes the turbibe to start 
rotating, and in a coil arranged around the axle, an alternating 
current is induced, whose frequency can be assigned to the 
rotation rate and hence, to m. It can be shown that this equip- 
ment - which functions satisfactorily for liquids (and large 
throughput rates - does no longer operate with the required 
precision in the case of small gas volumes, among other 
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reasons because of slippage and compressibility problems 
(see section 3.1). Furthermore, the calibration is possible 
only if an additional system is used (see point j), below). 

b) Pressure difference determination: 

The resolution of the Ap measurement range here demanded 

requires very small slit opening diameters (in the 0.1 mm 

range) and hence, is strongly dependent on the discharge 

coefficient 4», which is a function of combustion chamber 

pressure (or respectively, of the injection pressure: p + 

Ap. .). 

^inj 

Taking into account Mertel's work [101] and after some pre- 
liminary testing, this possibility was not further pursued. 

c) Perpendicular pipe with float element: 

A vertical glass tube with calibration marks is installed in 
the N 2 pipeline (flow direction: vertical, upwards), inside 
which a spherule can move, of a material selected as a func- 
tion of its specific gravity and the pipe flow’s dynamic 
pressure. Independently of reading problems, this method - 
often used with liquids - is not satisfactory, due to slip- 
page and environmental effects, for small gas quantities. 

d) Filling status measurement using a secondary medium (liquid): 
An alternative often used in technology is a liquid column 
connected to the gas storage container. The pressure remaining 
in the container can be estimated from the displacement of the 
liquid level, and hence, the quantity of gas used. Based on 
the required specification, the thermodynamic conditions (T-P 
problem) proved to be badly reproducible. An electronic indi- 
cation of liquid level is additionally impeded by the lack of 
adequate liquids (conducting, or optically evaluated, etc.). 

e) Measurement of the pressure increase upon discharge into a 
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large container: 

The possibility of determining m by recording the course of 
the pressure in a large chamber of known volume, as a func- 
tion of the valve setting of the rate of flow, had to be 
discarded, after intensive preliminary testing. This was due 
to the fact that because of the necessary valve opening times 
Joule-Thomson effects and environmental temperature influences 
could no longer be neglected. 

f) Direct optical measurement of the velocity of flow: /85 

This alternative also had to be ruled out, because the density 
differences in the N 2 working medium were not sufficient to 
render clearly visible changes in the flow velocities within 

the pipe. The admixture of particles that would positively 
change this characteristic had to be excluded because of the 
gas purity requirements. 

g) Direct inductive or capacitive measurement of the quantity 
discharged: 

Preliminary work was performed also on this concept, in con- 
nection with variation d), above. This alternative also had 
to be abandoned because of purity requirements (see also f)) 

h) Heated wire probes: 

After intensive collaboration with suppliers of such systems 
the integral measurements of the flow velocity and of the 
medium’s density (repurified N^) by heated wire anemometer 
procedures proved to be insufficiently precise for the small 
tube diameters here used (0 typicallv 2 mm) , as well as non- 
reproducible, because of capillary action and the influence 
of the probes on the flow itself. 

i) Laser Doppler velocity meter: 

Eventually, this laser-optical method, given appropriate 
development efforts on the systems available today, will permit 
direct measurements of the flow profile in narrow tubes and of 


87 


the nozzle discharge velocity distribution. Since currently 
the obtainable precision is still insufficient, the applica- 
tion of this method had to be disallowed, in the context of 
this study. 


j) Weighing system: 
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KEY I Gas consumption 2 Trapeze of equal surface area (approximated) 

3 Filling process 4 Nominal operation at fully formed ultrasound flow. The 
slight increase is due to gas cooling by the throttling effect 5 Emptying 
process 6 Opening the injector valve 7 Closing the injector valve 


Figure 44 Typical mass throughput course (time scale distorted) 


All methods described would require additional calibration 
methods in the event of their application, ultimately based 
on gravimetric methods. Hence the idea arises of using such 
methods directly, without the intermediate steps described 
above. These gravimetric methods can be divided into two 
main groups: 


- Continuous weighing -»■ m = f(t) 

- Integral weighing (before and after the flow process) 



"nach" = after 
"vor" = before 
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The tendency shown in Figure 44 (page 88) is decisive for 
the latter of these versions. 

The approximation of the 1^2 throughput is usually performed /87 
by application of the trapeze equation [68] 


r‘- r 

j ffl dt > A *1 


- t.) ♦ (t, - 


t,)\ 

J- », 


vor “ "nach 


(105) 




<%or ' %ach) • * 
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On the one hand, the trapeze method would be sufficiently 
precise (i.e., the filling and emptying processes could be 
neglected) if the increase in m due to the throttle effect at ’ 
very large differences ( t^ - t 2 ) did not occur. On the other 
hand, however, it is precisely for large values of ( t^ - t 2 ) 
that the throttle effect has an increasingly interfering 
effect, which forces keeping test times in the direction of 
smaller At values; because of the non-stationary gas consump- 
tion during switching on and off, this, in turn, worsens the 
precision . 

In case of continuous weighings of the N 2 consumption, these 
effects are of subordinate importance, since one can wait for 
stationary conditions before taking several m-values. This 
also makes it possible to evaluate drifts (caused by gas cool- 
ing) quantitatively. 


89 


3. 2. 3. 2 Selected weighing system - PC AAOO 

The main problem of this measurement method is ultimately the 
required resolution per weighing (better than 0.01 g), with 
the simultaneous utilization of a relatively heavy storage con- 
tainer (approximately 7-5 kg). The balance eventually decided 
on (Mettler, Zurich; PC 4A00), with a resolution of 0.001 g and 
a measurement range to A. A kg» is by far the best instrument on 
the market, for its range. Since no systems exists with c 
comparably high resolution beyond the measurement range of A. 4 
kg, the m-measurement problem becomes therfore displaced to the 
precompensation of the remaining 3.1 kg of the storage container 
(including the reduction station and the gas supply), which can 
not be balanced by design changes. Safety regulations, the test 
times required and the mechanical design overall do not allow 
for a container with a total mass substantially below 7.5 kg. In 
the required prebalancing, a decisive advantage of the basic 
balancing principle of the PC 4400 chosen, finds its expression: 
the mass to be determined presses on a vertically freely moving 
cylinder, kept in its position by the action of a strong magnetic 
field. An optical edge establishes this position to better than 
0.0001 mm and provide actual value information to a control 
circuit, such that it can change the coil current until the mag- 
netic force fully compensates the force due to the weight, and 
the cylindrical tod again occupies its original (unloaded). This 
means that the PC 4400 balance weighs without displacement. Hence 
the coil current is a direct measure of the corresponding mass. 

On this basis several prebalancing principles were examined 
theoretically and experimentally: 

- Counterweights with roller or knife-edge bearings 

In the end it was not possible to solve the mechanical bearing 
problems (close collaboration with the Chair for Precision In- 
struments, Prof. Unterberger [176]). 
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Upthrust body 

Immersion of hollow bodies into liquids or upthrust generation 
by means of gas-filled balloons had to be abandoned because 
of problems of reducibility due to environmental influences 
(mainly dependence of characteristic material values on 

Permanent magnet electromagnets /89 

These methods also had to be excluded because of their long- 
term behavior (decrease in field intensity) and dependence of 
the magnetic properties on following very extensive pre- 
liminary testing (in addition, the balance system was affected 
by the strong compensation fields required.) 

Spring system 

Since as we mentioned above, the PC AAOO weighs without 
displacement, spring compensation is appropriate 

^Feder “ ^ ^ ^ ("Feder” = spring) (107) 

The first preliminary tests already showed that small changes 
in the environmental temperature, near the spring - caused, 
for instance, by convection currents - changed the spring 
constant so sensitively, that reasonable measurements became 
impossible In order to solve this c = f(T^) problem, first a 
pendulum construction (see also clock pendulum) was implemented 
using three compression and three tension springs, in this case, 
made of the same material; its main characteristic is that it 
solves the temperature stability problem with respect to the 
relief spring [176]. Unfortunately the manufacturing problems 
in relation to a symmetrical force introduction in the six 
equally long spring proved to be too difficult. 

Hence we must forego a mechanical correction of the spring 
constant changes due to T^. The method chosen for use utilizes 
a cylindrical spring with extreme insulation thicknesses 
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surrounding the spring itself (without touching it!) and around 
the entire balance system (see also 3.2.2 and 3. A). In addition, 
this compensation spring does not act directly on the block 
but on a support thermally insulated from it. This prevents block 
shortenings (due to cooling effects of the gas in the reducing 
station) from affecting the sprin length and hence, the force /90 
of compensation. The final design is shown in Figure 32 and A?. 

In addition, the balance's power supply (energy loss heat) was 
provided with water cooling to maintain temperature gradients 
(convection) as low as possible. 



Figure 45 Balance printer control circuit diagram 

For a simpler output of the balance’s weight indications, these 
values are issued by a printer (Mettler GA AO). The printer 
system is activated by means of a quartz-stabilized control 
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designed by us. The circuit diagram for the printer control is 
shown in Figure A5, above. 

3. 2. 3. 3 Injector triggering /g^ 





Figure 46 Circuit diagram for the injection time control 

The control includes an adjustable time triggering (decadic sub- 

-4 IS 

division, from 10 to 10 sec and powers from 2 to 2 ). This 

triggering controls not only the printing of the instantaneous 

balance weight value, but also serves as time synchronization 

for the remaining measurement value acquisition and processing 

units (see also section 3.3). In addition, a time switch was 

developed to allow preselection of the test time and subsequently 

control the injection valve. In order to avoid feedback between 

the electromagnetically actuated valve and the electronics, the 

latter is electrically isolated from the valve’s actual control 

relais by an optic coupler. 

When the injector valve is activated, the p^ control of the 


93 






vacuum chamber is simultaneously switched off and full pumping 
power is released (see also section 3. 2. 2. 2). Because of the 
very precise test duration control it is also possible to 
control the continuous m-measurement here used by means of 
the integral method mentioned above. 



Figure 47 Mechanical assembly for throughput measurements (inside the 
balance's environmental chamber 
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The standardization and calibration of the entir m system 
was made possible using small steel beads (of known mass), 
which rolled into the receiving container, located in the 
isolated bottle block frame (see Figure A7, page 9^ ) » through 
a flexible tubing (i.e., the inverse of the procedure for 
discharge in the actual precision measurement test), at pre- 
cisely known time intervals. Standardization tests show that 
after solving the main problems (in relation to the spring 
constant c) the m-measurement construction operates very 
satisfactorily . 

3.2.4 Combustion chamber pressure determination 

In order to determine the combustion chamber pressure, an 
HMB-11 (Hottinger-Baldwin Co.) pressure sensor designed fol- 
lowing an inductive method is flange-mounted perpendicularly 
to the combustion chamber axis. It is connected to the chamber 
by a capillary (0 0.6 mm). 

The pressure to be measured acts on a thin membrane, which in 
turn hermetically seals off a reference volume at a pressure of 
1 bar. If this membrane is deformed due to a counterpressure, 
then it pushes a small piston attached to its center into an 
induction coil excited by a 100 kHz alternating current. Through 
a closed control circuit, the coil current is modified long 
enough (see section 3.2.3) - as a function of piston displace- 
ment - to return the membrane to its undeformed original posi- 
tion. After suitable electronic matching, the coil current 
allows a direct statement to be made regarding p^. In contrast 
to many of the customary systems, the advantage of this sensor 
is the measuring probe's constant dead volume. The precision 
stabilization of the P1 1 measurement tube supply voltage, the 
100 kHz carrier frequency generator, the membrane control cir- 
cuit, the measurement value amplification and its matching are 
all housed in a fully integrated instrument, the HBM-MC-IA. 
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In order to minimize interfering pick-up and power influences, 
the MC-IA electronics, housed in a metal cylinder, is attached 
directly to the PIT pressure measurement tube. The arrangement 
can be seen in Figure A3. It is thereby also possible to reduce 
the number of external incoming leads to five cables: two for 
the coarsely stabilized 12 V supply, one cable for screening 
or grounding, respectively and two for the measurement value. 

Since the pressure measuring device is mechanically connected 
to the combustion chamber and hence also to the thrust frame, 
these incoming leads must be fashioned of very thin wire (that 
is permissible, since only very low electric power is needed) 
and flexibly suspended. Outside of the test chamber, the 
measurement value must be amplified further, for a better 
application to the measurement value acquisition system. An /9A 
overview of the entire combustion chamber pressure detection 
system is shown in Figure A8, below, as a circuit diagram. We 



Figure 48 Measurement value amplification and adaptation of the combus- 
tion chamber pressure detection system 

should remark that the utilization of sensors based on the piezo- 
electric effect must be avoided, since they are not well suited 
to stationary pressures (as they occur here, essentially), in 
comparison to strong, dynamic pressure variations (the problem 
of long-term drift in charge amplifiers.) Careful calibration 
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tests have shown suitable agreement with the required speci- 
fications: reproducibility better than 0.3% of the measurement 
range of 1 to 25 bar absolute pressure. 

Due to the precise determination of the combustion chamber 
pressure it is possible to determine a 1-diroensional flow rate 
coefficient Cj., using a 1-dimensional mass throughput equation: 

, with the aid of the direct m-roeasurement method 
described in section 3.2.3 

D.t-d (J08) 

“c*” 

Typical results are shown in Section A. 

3.2.A.1 Determination of the p^-measurement tube installation /95 
error 

Departing from the simplified theory, the gas velocity inside 
the combustion chamber is not identically 0, but has a finite 
value. Because of this and determined by the technical boundary 
conditions due to the connection of a pressure measurement probe 
to a propulsion unit, the measured pressure deviates minutely 
from the stagnation pressure [12] 

J . ,2 (109) 

where = throat cross-section, 
c 

= installation cross-section 
M 

^ 6) ' 

B = angle between sensor axis and propulsion unit axis 

(here B = 0) 

With the boundary conditions applicable here, the pressure cor- 
rection is merely 
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"c • "c.Hess • 

and can therefore be neglected. 

In the course of the test execution it was confirmed that the 
initial assumption regarding the gas behaving adiabatically 
within the pipeline system and the combustion chamber, was 
correct. Because of the high heat capacity of the pipelines, 
made esssentially of copper, the cooling phenomena in the work- 
ing medium caused by throttle effects are compensated (provided 
care is taken not to extend the test duration to several minutes) 

3.2.5 Temperature measurements 


KD *>90 



Figure A9 Power supply and measurement value adapta- 
tion of the AD 590 temperature sensors 

In order to measure the N 2 gas temperature in the supply 
lines as well as in the combustion chamber of the cold gas pro- 
pulsion unit, and at several location of the test installation 
to control environmental conditions, temperature sensors Model 
AD 590 , of the Analog Devices company, are used. The decisive 
advantage of these models, in comparison to the customary thermo- 
elements consists in that both measurement amplifiers and supply 
voltage stabilizers - as well as a 0 °C reference temperature 
point - are integrated into the smallest colume (2 x A x 0.5 mm). 
The calibration and spreading followed the very simple circuit 
shown above. 
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The sensor proves to have a high resolution (better than 0.01 
°C) and an extraordinary long-term stability. The linearity 
error is less than 0.1% and dislocations due to contact 
problems, as they occur with thermoelements, have not been 
observed. Due to the low mass of the temperature-sensitive 
portion of the "sensor head", the response time (to 99% of the 
temperature value) is < 1 sec and hence certainly adequate to 
the application here contemplated. 


3.3 Data acquisition and data processing 
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All parameters essential to either measurement or test operation 

are recorded by high resolution analog recorders, for more rapid 

optical control (trouble shooting). Measurement variables that 

do not change during routine operation (such as lab’ lab’ 

p ) are rendered visible on meters. The variables necessary 
'^pump •' 

for performance determinations of the cold gas propulsion unit 
are in addition fed to the 10 channel ADC [analog to digital 
converter] developed at the Chair for Space Technology: it has 
an associated memory (capacity: 256 words 12 bits long = A decimal 
places). The circuit is shown in Figure 51 (page 101) (see also 
Roth, A. [127] and Neuwald, W. [111]). Following each measure- 
ment, the 10 data channels are printed out (including mean 
values) by means of a line printer. Through the use of modern 
microprocessors, the discretization errors of the ADC are at a 
level below' 0.1%, which can be disregarded. 

The block diagram below (page 100) provides an overview of data 
determination and output 

KEY to Figure 50 (page 100) I combustion chamber pressure 2 Measurement 
tube 3 Carrier frequency amplifier 4 Amplifier 5 Tracing 6 Relay 7 laser 
beam control 8 Motor drive 9 internal stabilization 10 Intensity control 
11 regulated laser 12 Thrust measurement 13 Measurement signal processing 
14 Active band filter 15 Low pass, active/passive 16 Thrust 17 Vacuum 
chamber pressure measurement and control 18 PID regulator 19 Vacuum chamber 
pressure 20 Magnetic valve 21 Needle valve 22 Vacuum pump suction line 
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Figure 50 


Block diagram of the data acquisition and 


processing method 
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(Key to Figure 50, continued) 

23 Mass throughput measurement 24 Balance (inductive) 25 Digitalization 
26 Printer 27 Evaluation 28 Temperature measurement 29 Calibration unit 
30 Temperature 31 Test duration control 32 to cold gas nozzle 33 Reducing 
station 34 Measurement value acquisition 35 lO-channel ADC 36 Analog 
control 37 Control unit 38 Memory 39 Computer (mean values) 



Figure 51 Block circuit diagram of the ADC and memory [27] 





A problem always to be solved anew in vacuum technology is that 
of electrical lead-ins through a vacuum chamber's walls. For the 
task here posed, 10 measurement channels and 15 supply lines 
(two of which carry 220 V AC) must be introduced into the chamber 
and remain isnulated from each other. Taking the cost problem 
into account, we opted for a newly developed 25-pole arrangement 
(manufacturing price: approximately DM 10.00). The savings per 
channel (in comparison to the customary 1-channel lead-ins such 
as those of the Leybold Heraeus company) is approximately 
DM 1,A00.00 


3. A Additional devices to reduce negative local influences /100 

The strong short-term variations in the laboratory line voltage 
(caused by the non-stationary power consumption of large external 
test facilities) force the utilization of the autonomous source 
of alternating power already described in section 3.2.2. 1 (see 
Figure 52). In order to reduce the temperature effects, espe- 
cially on the thrust and mass throughput measurements, the test 
arrangement is thermally isolated from the environment. To this 
end, a compelte "house” made from 10 cm thick styrofoam blocks, 
was built around the vacuum chamber and the laser device. In 
addition, temperature variations within the laboratory were 
reduced by insultaing the face of the windows and switching off 
the air conditioning unit, which operated much too irregularly. 

During "test-free" periods the "house", including all working 
openings, is kept open, thus reducing to a minimum the tempera- 
ture difference between the "house content" and the laboratory. 

Priot to the start of the test, the nozzle to be studied is 
built into the vacuum chamber and subsequently the thrust measure- 
ment chamber and the "house" are closed. The strongly insulating 
effect of the styrofoam insures that the effect of any small 
temperature differences between the vacuum chamber and the 
laboratory, during the test phase, can be ignored. In addition, 
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this coyering performs the functions of a wind screen, something 
absolutely essential because of the resolution required of the 

I 

thrust jne&suring device. 

' ' ■ . /1 01 



Figure 52 Independent 220 V AC power supply 



Influences due to stationary mechanical floor vibrations - such 
as those caused by a rotary sliding vane vacuum pump, for 
instance - can be canceled by using appropriately chosen mechan- 
ical dampers. Spontaneous strong disturbances such as those from 
traffic on the street, however, make it advantageous to perform 
any precision measurements during the night hours. 
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3.5 Construction of the cold gas propulsion unit 
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The measurements of the effect of e on the development of the 
boundary layer losses require the divergent nozzle length to 
be changed, without' affecting the remaining dimensions, 
especially in the throat area. The initially examined way of 
building the nozzle in slices, using thing thin foils of 
linearly increasing bore diameter, was impractical in its real- 
ization, since a multiplicity of very thin foild would have to 
be used, to avoid edge effects. In addition, accurate centering 
was not sufficiently possible. The alternative of boring the 
nozzle using precision drills or form cutters was closely 
examined and then discarded, following discussion with the 
appropriate companies. Neither the small bore dimensions nor the 
defined transition geometry in the nozzle throat area (backtaper) 
can be achieved in this manner. Arc erosive methods were also 
ruled out. While the craters produced thereby could subsequently 
be smoothed by grinding and polishing, satisfactory backtapering 
can not be obtained in this case, either. 

The galvanoplastic manufacturing method finally chosen uses a 
common core, fashioned in one piece, for both the convergent 
and the divergent portions of the nozzle. This core is machined 
"flying" on a small, precision turning machine (by a method for 
which a patent is currently being sought). Figure 5^ shows the 
main characteristics of this manufacturing step. 

The forming blades here used were sharpened by scanning a strongly 
enlarged former plate, which in turn was fashioned using a profile 
projector. (The manufacture was kindly undertaken by the Wendel 
Co. of Eichenau, near Munich.) 

Only the flying processing insures that the core does not break /10A 
at the narrowest cross-section (A^) - the narrowest value attained 
was 0.23 mm - since backrest clamping and secondary drive methods 
were eliminated due to exxcessively large misalignment). 

lOA 




Forming blad 




Figure 54 Schematic of the nozzle core's manufacture 

Following manufacture of the core, a 3 mm wall thickness is 
electrodeposited from a copper electrolyte bath. After removal 
from the plating bath, the external contour of the propulsion 
unit is fashioned in copper and subsequently the core is removed 
by gently twisting and carefully pulling apart at both ends. In 
order to make possible the simple mounting and demounting of the 
propulsion unit, each copper nozzle is fitted into an installa- 
tion flange and welded to it (for tightness). The individual 
steps in the nozzle manufacture are illustrated in Figure 55, 
page 106. Once a measurement series is completed, for a certain 
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Figure 55 Nozzle manufacturing steps (model embedded in casting resin) 

surface ratio e, the propulsion unit is converted - with the 
aid of a specially made holding device - to the next length or 
e step, by machining. 

Table ^ (page 10?) lists the main dimensions of the 12 propul- 
sion units examined. 



Figure 56 Original view of the 12 propulsion units examined 
(prior to first machining) 
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TABLE A MAIN GEOMETRIC DIMENSIONS OF THE 12 NOZZLES EXAMINED (NOZZLE NO. 1 WAS 
DAMAGED JUST PRIOR TO COMPLETING ITS MANUFACTURE) 


Nozzle 

Nr .: 

( mtn ) 

Bkt 

(mm) 

Bk ^ 

( mm ) 

(•) 

‘^200 

(mm) 

^200 

(mm) 

^200 

(-) 

°100 

(nrn) 

^00 

( mm ) 

^100 

(-) 

O 5 O 

( mm ) 

^50 

( mm ) 

^50 

(-) 

( mm ) 

'25 

( mm ) 

— 

'25 

(-) 

“to 

( mm ) 

^0 

( mm ) 

H 

2 

0,142 

— 

5,0 

2.0 

15,0 

3,89 

15,19 

187,89 

2,84 

13,22 

99.65 

1 . 9 ? 

11,47 

46,31 

1.41 

10,48 

24,59 

0,85 

9.48 

8,96 

3 

0,138 

5,0 

2.0 

20,0 

3,87 

13,67 

196,96 

2,80 

12,14 

102,73 

1,96 

10,96 

50,34 

1.37 

10,15 

24,81 

0,97 

9.44 

12.54 

4 

0,153 

3,663 

2.0 

10,0 

4.28 

20,78 

195,66 

3,09 

17,31 

101,64 

2.17 

14,68 

50,54 

1,51 

12,84 

24,53 0.93 

11.20 

9,30 

5 

0.176 

5.0 

2.0 

15,0 

4,92 

17,08 

194,06 

3,55 

14,48 

100,71 

2.54 

12,54 

51,88 

1.78 

11,09 

25,49 

t.n 

9,84 

9.95 

6 

0.166 

5.0 

2.0 

15,0 

4,68 

16,59 

198,07 

3,33 

14,05 

100,00 

2,39 

12,22 

51,68 

1,66 

10,84 

25,03 

1,03 

9.71 

9.68 

7 

0.165 

S.O 

2.0 

20.0 

4,59 

14,62 

192,39 

3,31 

12,85 

« 9,85 

2,40 

11,51 

52,89 

1,68 

10,49 

25,72 

1,00 

9,61 

9,18 

8 

0.206 

3,443 

2.0 

10,0 

5,79 

25,41 

196.58 

4,13 

20,70 

99,88 

2,93 

17,26 

50,58 

2,07 

14,78 

25,21 

1,29 

12,62 

9,82 

9 

0.215 

5,0 

2.0 

15,0 

5,87 

18,70 

186,36 

4,31 

15,77 

100,47 

3,07 

13,38 

51,11 

2,16 

11,61 

25,30 

1,34 

10,12 

9,73 

to 

0,223 

5.0 

2.0 

20.0 

6,18 

16,45 

191,75 

4,47 

14,11 

TOO . 00 

3,21 

12,30 

51,73 

2.25 

10,96 

25,36 

1,38 

9,79 

9,56 

n 

0.306 

3,092 

2.0 

10,0 

8,42 

30,94 

188,51 

6,09 

24,43 

98,54 

4,41 

19,45 

51,72 

3,08 

15,64 

25,22 

1,72 

11,85 

7,93 

12 

0.322 

5.0 

2.0 

15,0 

9,12 

24,40 

200,15 

6,51 

14,92 

101.87 

4.59 

15,77 

50,66 

3.22 

13,23 

24,95 

2,02 

11,03 

9,86 

13 

0,324 

5.0 



J .0 

20,0 

9,03 

20,24 

194,28 

6,53 

16,72 

101,47 

4,59 

14,03 

50,29 

3,25 

12,17 

25,13 

2.00 

10.51 

9.57 
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Microscopic examinations of the nozzle throat area before and 
after the test series showed that there was no substantial wear 
and hence, change in geometrical dimesnions. The evaluation was 
based on images from the MPI* Heidelberg electron microscope, 
provided by graduate engineer Peter Lell. Figure 57 shows the 
throat area, while Figure 58 represents a strongly magnified 
portion of the nozzle’s external contour. ^ 



Figure 57 Scanning electron 
microscope image of the throat 
area (Magnification, X 200) 




Figure 58 Partial 
view of the external 
contour 


* No expansion provided in the foreign original 


/106 


108 





k EXECUTION OF PERFORMANCE MEASUREMENTS 
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A.l Review of the test extension 

The twelve cold gas propulsion units were measured in five test 
series. After each of these series the nozzle's end cross-section 
was reduced by mechanical means (e-steps; 200, 100, 50, 25, 10) 
and measured by microscopy. Within a given test series, five 
combustion chamber pressures were applied four times each * for 
a total of 1200 individual tests. Few of the interferences with 
the measurement arrangement during a given test - caused mostly 
by eccentri loading of the m-balance of the N^-bottle frame - 
could be eliminated to complete the test. T'-e installation of a 
sufficient number and size of readily accessible work openings 
proved to be extraordinarily positive (see section 3.4). Due 
credit must be given the stabilization, in this sense. Any 
subsequent experimental efforts should include the following 
improvements, to increase precision and facilitate test execution: 

- fully airconditioned test space 

- separate data acquisition room 

- supply systems (mains operated power supply) outside of 
the test chamber 

- more stable building construction 

- low-traffic external environment 

See also section 3.4, in this regard. 

For a more ready control and handling of the large quantity of 
measurement values, in the subsequent calculations, we used 
the following nomenclature: 

* 20, 15, 10, 5, and 2 bar 
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Combu St ion 
chamber 

Nozzle number pressure E-step 

(coarse adjustment) 


Example 

2 
3 
A 

5 

6 


13 



S.10 

1 

0 1 X. 


02 

oi 200^ 

*^1 

05 

1 » 100 

2 

10 

2 « 50 

3 

15 

3 « 25 

4 

20 

4-10 

• 

• 


Sequential test number, 
for a given nozzle, 

'at the same pressure 
and at constant e 


4.2 Evaluation of measurement results by means of typical 
examples 
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The output from the measurement acquisition system arranges all 
experimental data according to the test number described above 
and makes them available to an evaluation program. Figure 78 
(see section 9-2) shows an input listing for such a test run. 
Here all measurement variables are examine for their order of 
magnitude and units, and corrected (for instance, the thrust 
measurement value is corrected for the Bourdon effect; the p^ 
value for p v, w ) • Next, especially tailored to the cor- 

responding test number, the theoretical performance prediction 
EGVP is started, such that in the end a data package arranged 
in parallel has been generated, with the theoretical and 
experimental performance parameters. In order to render the 
quantity of data thus created more tractable - by means of an 
interactive (i.e., terminal controlled) graphics program - 
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ISP-MES/ISP-0EL F-ME5/F-DEL 

!>,9 .« .* 1. 1 .0,0 .» .• t.J l.t 7,0 
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Pigure 61 Ratio of measured 
to predicted thrust, plotted 
as a function of e, with Pj. 
as' parameter 



) 

3) 


Pigure 63 Ratio of measured 
to predicted impulse, plotted 
as a function of e, with p^. 
as parameter 


KEY 1 Surface ratio e 2 Data package: result 31 3 Nozzle No. 4 from 


the data are plotted. Not only can the scale of resolution be 
arbitrarily changed, but it is also possible to perforn changes 
in the equalization curve, with respect to individual measure- 
ment points. Each of the 40 values per test run can be chosen 
to be X, y or a parameter point. Due to the interactive input 
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Eigure 62 Boundary layer loss 
as a function of the tbroat- 
Reynolds number, with e as 
parameter 



KEi I Boundary layer loss 2 Data package: result 31 3 Bozsle No. 4 from 

5 Throat Reynolds number 


from the graphics terminal, previous verification of the results 
can help in decisively reducing paper waste due erroneous output. 
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Figure 65 Mass throughput 
a function of with e as 
parameter 
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Figure 66 Measured over pre 
dieted mass throuput, with p 
as parameter * 


^5) cl_‘A^ct^E^fvER. EPSE 


les 


90S 


KEY 1 Data package: Result 31 2 

Figures 6l to 73 show typical 
evaluation. 


Nozzle No. 3 from 4 Pressure 
trends in the measurements 


5 E 

after 


In contrast to the 


corresponding diagrams of the 


m-ratio 





F-MES CNJ ^ MPKI-MES IK6/9J *10 




KEY 1 Data package: Result 31 2 Nozzle number 

^®mess^®del the^ > which shows good agreement between the measured 
and the predicted values of c^j over the entire range of e, we 
can see from Figures 61 and 63 (page 112) - which describe the 
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Y*1.4 

1.3 

1.2 


Pc'P. 

Figure 69 Flow separation ratio as a function of nozzle/environmental 

pressure 


thrust and impulse behavior, respectively - that there is a 
certain displacement in direction > 1.0 for the higher ranges 
of e. The explanation for this is as follows: the theory of 
performance prediction calculated here, at first, started from 
the premise (for lack of other information), that also for the /1 13 
case of smaller nozzles the so-called *'Summerfeld criterion” 
would be valid, with respect to flow separation during over- 
expansion. For the case p^<0.A p^, flow separation results 
and -*• impulse losses. Since in all 12 of the micronozzles 
examined the measured thrust at e > 100 was greater than for 
Summerfeld O.A (including all losses), we may infer that for 
the small propulsion units under discussion here, this factor 
is smaller than O.A. This means later separation of the flow 
from the wall and hence, a more widely usable e-range. Due to 
the very flexible evaluation program mentioned earlier, it is 
possible to change the Summerfeld factor until good agreement 
can be obtained. The results of that particular investigation 
show a value of up to 0.15. 
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Figure 71 Measured thrust as 
a function of predicted thrust 
with e as parameter 


C.C ,3 .♦ .6 .# 1.0 

SCHL'P-MESS IN N (1) 


KEY 1 Thrust measurement 2 Data package: Result 31 3 Nozzle No. 4 from 

The physical explanation of this discrpancy might be sought in 
that due to the greater influence of boundary layer formation 
in the case of extremely small propulsion unit dimensions (in 
contrast to full-scale propulsion units), the residual pressure 
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Figure 72 Throughput ratio as 
a function of combustion chamber 
pressure, with nozzle number as 
the parameter 



Figure 73 Specific impulse as 

a function of p ,, for 

°>,vac chbr 

the third estep 


KEY 1 Pressure 2 Data package: Result 31 3 e 4 from 5 Nozzle No. 

near the wall is higher than in purely potential flow. However, 
since separation occurs only when the wall pressure falls con- 
siderably below the environmental pressure (precisely for large 
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nozzles: factor O.A) outside of the nozzle, here this point is 
displaced backwards, in the direction of a larger surface ratio 
(see also [ 1 9 ,25 , A9, 90 , 95 , 108, 118, 122-12A , 139-1A3 , 1 A5 , 15A , 1573 ) . 


We find indications of a tendency to deviate froma factor of 
0.4 already in Kalt 180]. Figure 69 (page 116) shows a decrease 
in the separation factor to below 0.2 for increasing pressure 
ratio Pj-/p„ and decreasing isentropy exponent y. The expression 



P,, -6,2 



( 110 ) 


(corresponds to the solid 
line, in Figure 69) 


is a good approximation to the separation pressure ratios found. /114 

(/116) 


: result 31 
10 


0 0.5 1.0 1.5.10^ 

throat reynolds number 

Figure 74 Throat throughput coefficient c^ as a function of the throat 

Reynolds number 

The mass throughputs of the propulsion units and of the vacuum 

pumps, shown in these Figures, indicate a very precise agreement 

with the theory and provide a good control on the required pump 

constancy for different p , , values. This is the onlv 

“,vac chbr •' 
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way to verify a statement regarding stationary test conditions, 
as required here (and as we already assumed, at the start). 

The courses of the boundary layer losses as a function of throat 
Reynolds number, shown in Figure 62 (page 112), are in good 
agreement with the theory described in section 2.5» They also 
make it clear that in the case of micronozzles the bosses 

have a substantially larger effect (by a factor of up to 10) 
than, for instance, in booster propulsion units. 
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5 ERROR CALCULATION 
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Since in the determination of the theoretical/experimental 
characteristic magnitudes of a propulsion unit, necessary to 
a performance comparison, not all of them can be measured 
directly (for instance Igp» ^gff» t statement about 

the limits of error must involve Gauss' law of error propagation. 

Since the ODE functions are already quite complicated, is is 
obvious that their partial derivation - necessary for Gauss' 
method - will become rather cumbersome. Therefore, numerical 
results can no longer be obtained by quick estimates, but will 
require EDP [electronic ^ata processing] . On the basis of 
reading errors and the standard deviations of individual measure- 
ment variables, we obtain the following table of input quantities 
for the program (see page 122). 

The subsequent Figure 75 (page 123) shows a typical printout /1 18 
for a completetest . In addition to the Gaussian error, in each 
case also the maximum error is shown that results from linear 
superposition. We should add that the computation time (CPU 
[central processing unit] time required) for the error program 
is of the same order of magnitude as that required for the 
entire boundary layer program. 
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TABLE 5 LISTISG OF INDIVIDUAL ERRORS 


Balance 

^ 0.005 • 10"^ kg 

AP- 

5 1 X 

APc 

S 0.4 X 


^ 0,2 k 

AT. 

f 0.2 k 

Ay 

S 0.1 %. 

A^mol 

f 0.03 V., 

^ calib 

? 0.02 •/.. 

Aae 

< 0,1* 

Ac 

S 0,7 X 

A*-plpe 

< 0,02 X 

ATupthrust 

0.5 X will be compensated 

Am-,, 

«9S, in the thrust 

S 0,17 X 

frame 


A*^Bourdon 

S 1,0 Xwill be compensated(171 ) 

A^f(T) 

i 0.4 X 

AElectronics 

i 0,3 X 'for 30 sec test duration 

AF 

Z 

X 0.9 X 
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KEY 0 Illegible 1 Test number 2 Input data 3 Surface ratio 4 Environment 
pressure 5 Measured throughput 6 Measured thrust 7 Error 8 Function value 
9 Maximum error 10 Gaussian error 11 Mass factor 12 Thrust coefficient 
13 Effective velocity lA Characteristic velocity 15 Thrust 16 Specific 
impulse 

Figure 75 Error protocol printout for Test No. 9.2001 
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6 COMPARISON OF RESULTS FROM SECTIONS 2 AND A 
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The agreement between the predicted mass throughput coefficient 
Cp and that determined experimentally, is very satisfactory 
(see Figures 72 and ?A ) . The predicted specific propulsion unit 
impulse data (especially I^p ^el^ have been clearly confirmed 
by measurements, in the range of e = 10 to 100 and p^ = 2-10 
bar. Beyond e = 100 and p^ = 10 bar, some uncertainty should be 
assumed in the overall statement, because Pg<<Poo snd due to 
the associated flow separation effects. We expect an improvement 
in future experiments, using vacuum pumps of higher suction 
capacity and developing a new method for the direct determina- 
tion of the separation site on the nozzle wall, and of the 
separation pressure, P___,. From the point of view of engineering 

O w ^ 

manipulation, the model used for boundary layer calculations may 
be considered satisfactory, since it allow sufficient prediction 
precision over a wide range of propulsion unit sizes and opera- 
ting conditions. Taking into consideration the simplified 
equation discussed in section 2.2.10, a fast answer in regard 
to can be obtained in subsequent applications. Further 

improvements seem possible, in the future, in the measurement 
system; in connection with improved computer programs it should 
become possible to further improve the reliability of the 
statements obtained using 2.2.10. The results obtained in the 
present case, for stationary operation behavior, can also be 
extended to pulse operation (as is often the case in attitude 
control tasks). The digital acquisition of measurement values 
and their subsequent processing and comparison with theoretical 
predictions represents a considerable simplification of efforts 
and makes it possible to eliminate large sources of error often 
encountered in manula calculations. 
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7 SUMMARY AND OUTLOOK 
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It is possible, using the nozzle manufacturing method described, 
to produce not only conical nozzles with throar diameters in 
the 0.1 mm diameter range, but also optimized "bell nozzles", 
with specific impulses up to 2% higher. It can be shown that 
their application in future satellite systems (especially if 
the lifetime is > 5 years) can result in the saving of several 
kg of fuel (for a 1 ton satellite mass, for instance), which in 
turn increases the payload capacity. Since as of this date it 
is hardly possible to further reduce the structural factors, 
and a limit is going to be reached in the near future also for 
specific electronic masses, a decrease in the orbit and attitude 
adjustment fuel requirements is absolutely necessary. 

The non-contact thrust measurement procedure here developed can 
be utilized also in other areas, in which the precise measurement 
of not directly accessible test objects is necessary. 

The use of non-displacement weighing systems with simultaneous 
prebalancing of the storage container is an area usable in measure- 
ment technology (resolution: nearly six orders of magnitude of 
contnuous mechanical processes), which is urgently needed to 
continue development of precise throughput measurement instruments 
and their calibration. By using the impulse-specific characteristic 
variables determined in the precision measurements, it is possible 
no only to more precisely adapt the extensive explicit computer 
programs, but also those more easily handled in the rapid design 
of propulsion units, providing "rule of thumb" but usable equa- 
tions. In the continuation of these performance measurements 
using larger vacuum pumps (higher suction capacity), these results 
can be extended to large (>150) nozzle surface ratios. 

The results of this effort can be used not only in space travel 
per se, but conceivably also, in the development of gas dynamics 
laser systems or the U-235 enrichment by nozzle-separation methods. 
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9 APPENDICES 
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9.1 Design of the cantilever beam 



KEY 1 Bending carrier 2 Mirror 3 Laser instrument A Array diode 

Figure 76 Arrangement of the bending carrier with mirror, laser instrument 
and array diode 

Specified quantities: 


G 

L 

max 

S 


a 


b 

max 

AS 


15 N (weight of nozzle and harness) 
250 mm 
0.5 m 

1 N (thrust) 

1 0 ym 
1 . 5 mm 
10“^ N 


For the bending line w(x) we have, for small deformations with 
a line load q due to the carrier’s own weight: 


E I w 


IV 


q 
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Following quadrouple integration 


E I w 


1 

* - wn 
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With the boundary conditions 


w{0) • w'iO) • W(L) - 0 


the integration constants become 


C3 - 0 
C 2 ■ 0 

- L F ♦ 2 

We then have, for the bending line and its inclination ■ 

*f(x) ■ i *^*^ ■*^ 7 (L F ♦ 2 

w'(x) ■ r * 7 ♦ X (L F ♦ ^ q)^ 

At the location of the mirror, x = L and we have 


w'a) 


ITT 


(here we assume q = 0 and F = G + S) 


The angular change of the reflected laser beam will be twice the 
change of the slope in the bending line, due to the laws of 
reflection; 


a 


[■' 


(L) 


G-^S 



/1A0 


The path b traversed by the laser beam’s point of impact on the 
array diode is 


, b « tan a 
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Since b is very small compared to L and therefore a is very 
small, tan a = o and hence, 


b • -pp • S 

The carrier must satisfy the following conditions: 

1 ) The deformation must be large enough to attain the desired 
precision 


‘■s *■ 

2 < g I '" • AS 


2) The deformation must be smaller than the size of the array 
diode allows 

L 

b > S 
max E I max 


3) No plastic deformation of the tube under thrust or its own 
weight, respectively 


o < o 


zul 


From 1) and 2): 


^max 


> *■ > a 


with 


o ■ . f_k 

max Wjj ^Itrelsrohp * ”*'"0 


(Kreisrohr = circular tube) 


/Ul 
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or 




> 200 mm 
• 10"^ N 
■ 0,5 ffl 

. N 

mm* 


as well as 

{Anticipating the result) 
("Stahl" = steel) 


we obtain: I 


100 mm 


A. 


if a 1 mm wall thickness is chosen, then 



or 


0 < C 

max 


r2/3 


From this we realize that L must be chosen as large as possible /1A2 
(within the vacuum chamber) . 


For the certainty against plastic deformation we have 

o. 


^plast ^ o. 


max 


With the design magnitudes listed above, becomes = 2.7, 

i.e., no weight compensation is necessary I 
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Of the materials examined, steel had the best qualities, in ad- 
dition to its preference because of ready availability and simple 
processing (can be welded). 

The nozzle holding flange welded to the free end of the cantilever 
beam (tube), including p tube (see Figure 77 , page 150) remains 

w 

clearly below 15 N, in weight, in the final design; hence the 
design calculations remain valid (with respect to S , ) . In 

P J. 3 S L 

addition the resolution of the array diode system could be con- 
siderably increased (see section 3.2.3). 
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P^-Druck- 

messbohrung 


SerrtffA-ff 
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\0 Bohrungen dfenen der Arrttl&rvng 
des KragbaUers wShrend des Trieb- 
werkaustausches 


KEY Detail X ^ preaaure measurement boring 3^ Section A Thruster 5 p meas. tube receptacle 
^ Part 3 View Y 2. ^2 Inflow line Part 2 9 ^ [Illeg.] JjO Boreholes serve to arrest the cantilevered 
beam during the propulsion unit exchanges 

Figure 77 Thruster receiving flange 






9.2 Flow diagram for computer program EGVP 
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KEY 1 Thermodynamic input 2 Geometry input 3 Control input 4 Calcula- 
tion of combustion chamber values 5 Boundary layer thickness calculations 
6 Core flow calculations 7 Geometry calculations 8 yes 9 Mach number cor- 
rection 10 No 11 End value calculation 12 Loss calculations 


Figure 78 Computer program flow diagram 
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KEY 1 Test No. 2 Date 3 Bottle 4 Combustion chamber 5 Vacuum 6 Room 7 Balance 
8 N 2 chamber. 9 Thrust 10 Remarks 

Figure 79 Example of test protocol (5.1001) 
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9 . A Cost analysis (180 ♦ 188) 


Balance and printer (Mettler) 


DM 7,700 


Laser lOlP (Laser Optronic) 

Components for measurement data acquisition 


1 ,A00 

600 (approx) 
A50 

PIN-LSC array diode system (incl. evaluation syst.) 1,800 


measurement tube and amplifier 


Neutral wedge 
Micrometer drive motor 
Vacuum minimotor 

Temperature measurement sensors ( 6 ) 
Molding knife for nozzle manufacture 
Nozzle manufacture ( galvanoplastics ) 
N 2 - 2 liter bottle ( 2 ) 

1 transfer tube (special manufacture) 
Small parts (approx.) 

Working gas (approx.) 

Documentation (approx.) 


190 
80 
140 
180 
1 ,500 
280 
300 
210 
500 
600 
700 


Total (approx. ) 


DM 16,630* 


* Of this, approximately DM 11,500 are outside financing 
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